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METHOD AND SYSTEM FOR PROCESSING 
POSITIONING SIGNALS IN A STAND-ALONE MODE 



RELATED APPLICATIONS 

This application is a divisional of U.S. Patent 
Application Serial No. 09/908,011, filed July 18, 2001 by 
Louis Henry Martin Jandrell and ENTITLED "METHOD AND 
5 SYSTEM FOR PROCESSING POSITIONING SIGNALS IN A STAND- 
ALONE MODE . " 

This application is related to copending U.S. Patent 
Application Serial No. 09/909,274, filed July 18, 2001, 
entitled "METHOD AND SYSTEM FOR PROCESSING POSITIONING 

10 SIGNALS WITH MATCHING ASSISTANCE" (Attorney's Docket 
068796.0102); copending U.S. Patent Application Serial 
No. 09/908,403, filed July 18, 2001, entitled "METHOD AND 
SYSTEM FOR PROCESSING POSITIONING SIGNALS IN A GEOMETRIC 
MODE," now U.S. Patent No. US 6,515,620 Bl issued 

15 February 4, 2003 (Attorney's Docket 068796.0104); and 
copending U.S. Patent Application Serial No. 09/908,199, 
filed July 18, 2001, entitled "METHOD AND SYSTEM FOR 
DETERMINING CARRIER FREQUENCY OFFSETS FOR POSITIONING 
SIGNALS," now U.S. Patent No. US 6,529,160 B2 issued 

20 March 4, 2003 (Attorney's Docket 068796.0105). 

TECHNICAL FIELD OF THE INVENTION 

This invention relates generally to the field of 
geolocation and more particularly to a method and system 
25 for processing positioning signals in a stand-alone mode. 
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BACKGROUND OF THE INVENTION 

In the application of global positioning system 
(GPS) technology to the geolocation of wireless devices, 
a typical geolocation function utilizes a course 
5 acquisition (C/A) code, or Gold code, which is received 
repeatedly from GPS satellites, in order to determine 
position. In addition to the repeated Gold code 
sequence, the Gold code comprises satellite message data 
that is modulated on top of the Gold code signal by 

10 inverting the phase according to the message data. 

The spectral density level of a signal from a GPS 
satellite received at a conventional GPS receiver with a 
direct line-of -sight to the satellite is significantly 
less than the thermal noise level of the conventional GPS 

15 receiver. When satellite signals are received at 

wireless devices being operated with obstructed views of 
the sky and thus obstructed line-of -sight , the satellite 
signals are weakened even further. Consequently, the 
obstructed signal levels from satellites are generally 

20 well below the threshold at which receivers may receive 
reliable message data signals from the satellites. 

Recent solutions to the problem of receiving 
weakened positioning signals provide for partitioning the 
geolocation processing functions such that some of these 

25 functions are performed at the unknown location which is 
to be determined and other functions are performed at a 
location with an unobstructed view of the signal source. 

For example, one of these methods, using a satellite 
signal source, provides for measuring all the satellite 

3 0 signal parameters, including the Doppler shift for each 
satellite signal, at unobstructed receivers located near 
the unknown location. The unobstructed receivers then 
send pertinent data to the unknown location to allow 
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p.eudorange estimation to be completed at the unknown 
location. 

Disadvantages associated with this example include a 
relatively expensive requirement of integration of the 
network of such receivers with the wireless carrier 
network that provides a link between the unknown location 
and the unobstructed receivers. In addition, this 
solution restricts the joint operation of the unknown 
location and the unobstructed receivers in performing the 
geolocation processing functions to only those carrier 
networks that are so integrated. 
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SUMMARY OF THE INVENTION 

The present invention provides an improved method 
and system for processing positioning signals in a stand- 
alone mode. This invention substantially eliminates or 
5 reduces disadvantages and problems associated with 
previous systems and methods. In a particular 

embodiment, the time to estimate a pseudorange from 
received positioning signals is minimized, while the 
processing gain is maximized to facilitate rapid 

10 detection of positioning signals while minimizing the 
consumption of energy. 

In accordance with one embodiment of the present 
invention, a method for processing positioning signals in 
a ranging receiver in a stand-alone mode is provided. 

15 The method includes collecting pseudorange samples from 
positioning signals received at the ranging receiver from 
a" plurality of satellites. The pseudorange samples 
comprise message data modulation. Each satellite has an 
associated Gold code. A previously determined carrier 

20 frequency offset (CFO) is selected from a plurality of 
directly extracted CFOs . The pseudorange samples are 
compensated for the selected CFO. The message data 
modulation is removed from the pseudorange samples . The 
pseudorange samples are stacked for each satellite. The 

25 Gold code associated with each satellite is correlated to 
generate a pseudorange time sequence for the satellite. 
A determination is made regarding whether an adequate 
correlation peak exists in each pseudorange time 
sequence. A pseudorange is determined for the ranging 

3 0 receiver based on the correlation peaks when an adequate 
correlation peak exists in each pseudorange time 
sequence . 
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Technical advantages of one or more embodiments of 
the present invention include providing an improved 
method for processing positioning signals in a stand- 
alone mode. In particular, the time required to obtain 
5 pseudorange estimates from weakened positioning signals 
is minimized and the processing gain available in any 
sample segment used to obtain pseudorange information is 
maximized. In addition, by using a direct extraction 
method to quickly determine a CFO, resolving the identity 

10 of satellites whose signals are being received, and 
compensating the received signal by the amount of the 
CFO, the pseudorange estimation may be accomplished 
without performing a time-consuming search through the 
Doppler frequencies and satellite codes. 

15 Other technical advantages of one or more 

embodiments of the present invention include a 
geolocation processor that may supply, to a receiver, 
current fragments of message data that were transmitted 
during the time the receiver -was collecting signal 

2 0 samples. This is possible because the geolocation 

processor and the receiver are able to exchange time-of- 
day synchronizing information. As a result, the 

modulation for the message data may be removed from the 
samples being processed by the receiver. Accordingly, 
25 the available processing gain that can be achieved from a 
signal averaging process is maximized. 

Due to the improvement in processing gain, technical 
advantages of one or more embodiments of the present 
invention also include an ability either to obtain 

3 0 greater sensitivity in a given amount of signal 

processing time or to significantly reduce the amount of 
processing time to reach a particular level of 
sensitivity. In this regard, a dynamic process is used 
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that acquires only the amount of signal necessary to 
achieve reliable detection. This minimizes both the 
processing time and the amount of intermediate -result 
memory storage required during signal processing. 
5 Yet another technical advantage of one or more 

embodiments of the present invention includes the removal 
of a requirement to integrate a geolocation processing 
system with the wireless carrier network that provides a 
communication link between the ranging receiver and an 

10 assisting function. Thus, the need for, and the costs 
associated with, a wireless carrier or plain old 
telephone system network are eliminated. For example, 
support from a carrier's network is not required to 
determine a coarse position estimation. In addition, a 

15 support network of nearby reference receivers, or its 
equivalent, is not required to provide Doppler-shif t 
compensation or Doppler-shif t search assistance. An; end- 
to-end frequency-calibration embodiment in the wireless 
communication link between the receiver and the 

20 geolocation processor to allow making use of the Doppler- 
shif t is also not required. 

In addition, technical advantages of one or more 
embodiments of the present invention include minimized 
memory use, power consumption and network loading, 

25 increased sensitivity, decreased amount of time for a 
location estimate, and greater privacy for customers. 

From the following figures, description, and claims, 
other technical advantages of the present invention will 
be readily apparent to one skilled in the art. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

For a more complete understanding of the present 
invention and its advantages, reference is now made to 
the following description taken in conjunction with the 
5 accompanying drawings, wherein like reference numerals 
represent like parts, in which: 

FIGURE 1 is a block diagram illustrating a 
geolocation system operable to process positioning 
signals in accordance with one embodiment of the present 
10 invention; 

FIGURE 2 is a block diagram illustrating the ranging 
receiver of FIGURE 1 in accordance with one embodiment of 
the present invention; 

FIGURE 3 is a block diagram illustrating the ranging 
15 receiver of FIGURE 1 in accordance with another 
embodiment of the present invention; 

FIGURE 4 is a flow diagram illustrating a method for 
obtaining a calibration location estimate for the ranging 
receiver of FIGURES 1, 2 or 3 in accordance with one 
20 embodiment of the present invention; 

FIGURE 5 is a flow diagram illustrating a method for 
obtaining a location estimate based on a non- initial 
trigger for the ranging receiver of FIGURES 1, 2 or 3 in 
accordance with one embodiment of the present invention; 
25 FIGURES 6A and 6B are flow diagrams illustrating a 

direct extraction method for determining carrier 
frequency offset in the methods of FIGURES 4 and 5 in 
accordance with one embodiment of the present invention; 

FIGURE 7 is a flow diagram illustrating a method for 
30 determining pseudorange in the ranging receiver of 
FIGURES 1, 2 or 3 in a stand-alone mode in accordance 
with one embodiment of the present invention; 
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FIGURES 8A and 8B are flow diagrams illustrating a 
method for aligning message data to determine a 
pseudorange in the ranging receiver of FIGURES 1, 2 or 3 
in a stand-alone mode in accordance with one embodiment 
5 of the present invention; 

FIGURES 9A and 9B are block diagrams illustrating a 
system for generating sub- stacks and stacks to determine 
coarse satellite message data alignment in the system of 
FIGURE 1 in accordance with one embodiment of the present 
10 invention; 

FIGURE 10 is a block diagram illustrating a system 
for refining satellite message data alignment in the 
system of FIGURE 1 in accordance with one embodiment of 
the present invention; 
15 FIGURE 11 is a flow diagram illustrating a method 

for aligning message data in the ranging receiver of 
FIGURES 1, 2 or 3 in accordance with one embodiment of 
the present invention; 

FIGURE 12 is a flow diagram illustrating a method 
20 for predicting satellite message data received in 
positioning signals at the ranging receiver of FIGURES 1, 
2 or 3 in accordance with one embodiment of the present 
invention; 

FIGURES 13A, 13B and 13C are flow diagrams 
25 illustrating a method for determining a pseudorange in 
the ranging receiver of FIGURES 1, 2 or 3 in a geometric 
mode in accordance with one embodiment of the present 
invention; 

FIGURE 14 is a flow diagram illustrating a method 
3 0 for determining geolocation in the ranging receiver of 
FIGURES 1, 2 or 3 with matching assistance in accordance 
with one embodiment of the present invention; and 
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FIGURE 15 is a flow diagram illustrating a method 
for determining geolocation in the geolocation processor 
of FIGURE 1 by providing matching assistance .n 
accordance with one e^odiment of the present invention. 



5 
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DETAILED DESCRIPTION OF THE INVENTION 

FIGURE 1 is a block diagram illustrating a remote 
geolocation system 10 in accordance with one embodiment 
of the present invention. The system 10 may comprise a 
5 global positioning system (GPS) , a global navigation 
satellite system (GLONASS) , or other suitable navigation 
system . 

In the case of GPS, satellites in a navigation 
constellation all transmit on the same 1575.42MHz 

10 frequency and use a code division multiple access (CDMA) 
technique to distinguish one satellite's signal from the 
others. In the case of the GPS constellation, a 1023- 
chip Gold code is biphase -modulated on the sinusoidal 
carrier. The Gold codes, unique to each satellite, are 

15 used for the CDMA and ranging functions. The chipping 
rate in this case is 1023 kHz, so that the transmission 
of the Gold Code repeats every millisecond. One cycle 
through the Gold Code is called a "code epoch." A set of 
samples whose duration is equal to the duration of a code 

20 epoch is referred to herein as an epoch. Navigation 
message data is modulated onto this Gold code modulated 
carrier signal by either inverting or not inverting the 
phase of the Gold code in sympathy with the values of the 
navigation message data bits. Each transmitted 

25 navigation message data bit coincides with 20 Gold code 
epochs, or 20 milliseconds. 

The system 10 comprises a plurality of positioning 
signal transmitters 12 and at least one ranging receiver 
14. The positioning signal transmitters 12 are operable 

30 to transmit positioning signals, and each ranging 
receiver 14 is operable to receive the positioning 
signals and to perform geolocation processing functions 
in order to determine a pseudorange estimate for the 
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ranging receiver 14 based on the received signals. As 
used herein, "each" means every one of at least a subset 
of the identified items. In the illustrated embodiment, 
the positioning signal transmitters 12 comprise 
5 satellites 12. Thus, as used herein, "satellites" means 
conventional satellites, pseudolites, or other suitable 
positioning signal transmitters. 

As described in more detail below, a location 
estimate, or geolocation, may be determined for the 

10 ranging receiver 14 based in part on the pseudorange 
estimate. According to one embodiment, a location 
estimate provides a location for the ranging receiver 14 
that is accurate to within 10 to 50 meters or other 
suitable distance. 

15 A location estimate may be determined when the 

ranging receiver 14 is receiving weakened positioning 
signals, such as signals weakened by obstructions. 
Weakened positioning signals comprise signals that are 
received from the satellites 12 at levels that result in 

2 0 insufficient information being received by the ranging 
receiver 14, in the absence of other satellite signals or 
other sources of information, to determine a location 
estimate but sufficient information being received to 
estimate a pseudorange for the ranging receiver 14 . 

25 According to one embodiment, weakened positioning signals 
comprise signals compromised up to 16dB below "clear 
view", or unobstructed, signal levels. In an alternative 
embodiment, weakened positioning signals comprise signals 
compromised up to 30dB below "clear view" signal levels. 

30 The system 10 also comprises a wireless device 16 

and a service center 18 which may be coupled to each 
other through a wireless carrier network 20 and the plain 
old telephone system {POTS) network 22. In addition, the 
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system 10 comprises a geolocation processor 24, which may- 
be coupled to the service center 18 either directly or 
indirectly, such as through a network or other suitable 
form of indirect communication and may be coupled through 
5 the service center 18 to the wireless device 16. 
Although the illustrated embodiment includes one service 
center 18 and one geolocation processor 24, it will be 
understood that the system 10 may comprise any suitable 
number of the components 18 and 24 without departing from 

10 the scope of the present invention and that functionality 
of the system 10 may be otherwise suitably distributed or 
combined. It will be further understood that a single 
geolocation processor 24 may be coupled to a plurality of 
service centers 18 without departing from the scope of 

15 the present invention. 

The wireless device 16 may be coupled to the ranging 
receiver 14, which may comprise push buttons 30, or other 
suitable interfaces, for user operation such as to 
initiate a location estimating process for the ranging 

20 receiver 14. The wireless device 16 is operable to 
provide a communication link between the ranging receiver 
14 and a service center 18 via the wireless carrier 
network 20 and the POTS network 22 and to provide a 
communication link between the ranging receiver 14 and a 

25 geolocation processor 24 via the service center 18. 

The wireless device 16 may comprise a cellular 
telephone, a wireless Personal Digital Assistant, a two- 
way or response pager, a private or commercial vehicle 
tracking system, an "On-Star" type motorist service 

30 network, a private or commercial wireless data network 
(or a device in such a network) , or any other suitable 
device or network capable of providing a data- link 
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between the ranging receiver 14 and the geolocation 
processor 24 . 

The wireless device 16 may be coupled to the ranging 
receiver 14 in any suitable fashion. For example, the 
5 wireless device 16 and the ranging receiver 14 may be 
coupled together by a cable 32. Alternatively, the 
ranging receiver 14 may be located near the wireless 
device 16 and linked to the wireless device 16 by a low- 
powered short-range wireless link, such as those provided 

10 by the Blue-Tooth protocols. The ranging receiver 14 may 
also be coupled to the body of the wireless device 16 as 
a "stick-on" attachment, or the ranging receiver 14 may 
be built into a receptacle such that the wireless device 
16 is in electrical contact with the ranging receiver 14 

15 when the wireless device 16 is inserted into the 
receptacle . 

In another alternative, the ranging receiver 14 may 
be integrated with the wireless device 16 . For this 
embodiment, the wireless device 16 is operable to process 

20 communication at a minimum of two frequencies, one 
frequency for receiving signals to be processed by the 
wireless device 16 independently of the functionality of 
the ranging receiver 14, such as cellular signals for a 
wireless device 16 comprising a cellular telephone, and a 

25 second frequency for receiving positioning signals. The 
wireless device 16 is also operable to digitize samples 
at a rate corresponding to the rate at which the ranging 
receiver 14 may process samples in accordance with the 
present invention. For example, the wireless device 16 

30 may be operable to digitize samples at a rate of at least 
2.048 Megasamples/second at a quantizer resolution of a 
minimum of two bits. It will be understood that other 
suitable rates may be used, as described in more detail 
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below in connection with FIGURE 2. In addition, the 
wireless device 16 may include additional memory, as 
compared to a wireless device 16 not coupled to a ranging 
receiver 14, in order to store samples for processing 
5 during pseudorange estimation. 

The service center 18 may comprise a Public Safety 
Access Point (PSAP) , a Vehicle Dispatch Application, or 
other suitable entity capable of facilitating geolocation 
processing for the ranging receiver 14 by coupling the 
10 ranging receiver 14 to the geolocation processor 24 . The 
service center 18 is also operable to manage the 
detection, routing and switching of voice traffic that 
may be associated with the geolocation determination 
process . 

.15 The service center 18 includes a rack 4 5 of modems 

40, a switch 42, and at least one telephone handset 44. 
Each modem 40 may include a Location Services Protocol 
modem or other suitable modem and is operable to provide 
a communication link between the service center 18 and 

20 the ranging receiver 14 associated with the wireless 
device 16 through the POTS network 22 and the wireless 
carrier network 20. The modem 40, or other communication 
link such as an Ethernet, may also provide a link between 
the service center 18 and the geolocation processor 24 . 

25 Each modem 4 0 is also operable to route voice 
communication to the switch 42 and data communication to 
the geolocation processor 24. The switch 42 is operable 
to route voice communication from a modem 40 to a 
telephone handset 44 and from a telephone handset 44 to a 

30 modem 40. The telephone handset 44 allows an operator at 
the service center 18 to interact with a user of the 
wireless device 16 through voice communication provided 
through a modem 40. 

DALOl: 756059.1 



ATTORNEY'S DOCKET 
068796 . 0103 



15 



PATENT APPLICATION 



The geolocation processor 24 is operable to identify 
the satellites 12 by matching the differences between 
pairs of carrier frequency offsets measured at the 
ranging receiver 14 with the differences between pairs of 
5 Doppler frequency shift estimates made of satellites in 
view at particular locations on the Earth's surface, to 
accurately calibrate a real-time clock in the ranging 
receiver 14, to track satellite message data from the 
satellites 12 and bit transition times in the satellite 

10 message data, to provide fragments of the satellite 
message data from the identified satellites 12 to the 
ranging receiver 14, to convert a pseudorange estimate 
from the ranging receiver 14 into a geolocation of the 
ranging receiver 14, and to provide the geolocation to a 

15 requester such as the ranging receiver 14 or PSAP 
operator. 

The geolocation processor 24 comprises a workstation 
46 and a receiver 48. The workstation 46 is operable to 
perform geolocation processing functions. The 

20 workstation 46 is also operable to obtain current 
ephemeris, time, and other suitable satellite message 
data from the receiver 48. Ephemeris data may comprise 
trajectory data received from one or more satellites 12 
that allow positions of satellites in the satellite 

25 constellation to be predicted. The receiver 48 may 
comprise a GPS receiver or other suitable receiver 
operable to receive positioning signals from the 
satellites 12. 

The wireless carrier network 2 0 may be operable to 

30 coordinate location management, authentication, service 
management, subscriber management, and any other suitable 
functions for a plurality of wireless devices 16. The 
wireless network 20 may comprise a public land mobile 
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network or other suitable network operable to provide 
communication between the wireless devices 16 and the 
POTS network 22 . 

The POTS network 22 may be coupled to the wireless 
5 carrier network 20 and to the service center 18 through 
communication links 50. The communication links 50 may 
be any type of communication link capable of supporting 
data transfer. In one embodiment, each of the 

communication links 5 0 may comprise, alone or in 

10 combination. Integrated Services Digital Network, 
Asymmetric Digital Subscriber Line, Tl or T3 
communication lines, hard-wire lines, telephone lines or 
wireless communication. It will be understood that the 
communication links 50 may comprise other suitable types 

15 of data communication links. The communication links 50 
may also connect to a plurality of intermediate servers 
between the POTS network 22 and both the wireless carrier 
network 2 0 and the service center 18. 

The system 10 may also comprise one or more remote 

20 receivers 60, service requesters 62 and/or application 
servers 64 which may be coupled to the geolocation 
processor 24 through a data network 66 . The remote 
receiver 6 0 may comprise a GPS receiver or other suitable 
receiver operable to receive positioning signals from the 

25 satellites 12. Data derived from these positioning 
signals may then be provided to the geolocation processor 
24 through the data network 66. Thus, the geolocation 
processor 24 may receive positioning signals from either 
receiver 48 or receiver 60. According to one embodiment, 

30 the system 10 comprises three or four remote receivers 60 
spaced around the Earth such that, at any given time, the 
three remote receivers 60 may receive current ephemeris 
data from each satellite in the constellation. 
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The service requester 62 may comprise a location 
information requester that is operable to generate a 
request for a location estimate for a particular ranging 
receiver 14 and to provide the request to the application 
server 64 through the data network 66. The application 
server 64 may comprise a server for location-based 
applications or other suitable server that is operable to 
receive the request generated by the service requester 62 
and to process the request. 

In processing the request, the application server 64 
is operable to provide the request to the. geolocation 
processor 24 through the data network 66, prompting the 
geolocation processor 24 to initiate a location 
estimating process for the ranging receiver 14 identified 
by the service requester 62 in the request. In addition, 
the application server 64 is operable to receive a 
resulting location estimate from the geolocation 
processor 24 through the data network 66 and to provide 
the location estimate to the service requester 62 through 
the data network 66. 

The data network 6 6 may comprise a public switched 
telephone network, an integrated services digital 
network, a local area network, a wide area network, a 
global computer network, such as the Internet or other 
dedicated switched network, or other communication system 
or combination of communication systems at one or more 
locations. According to one embodiment, any of the 
networks 20, 22 and 66 may share a same physical 
implementation with one or both of the other networks 20, 
22 and 66. 

Because the Doppler frequency shifts are 
continuously changing due to the movement of the Earth 
and of the satellites 12, the workstation 46 may be 
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Operable to regularly re-compute the pattern of Doppler 
frequency shifts of satellites 12 based on the changing 
satellite positions for use in satellite identification 
by Doppler frequency shift differences, as described in 
5 more detail below. The workstation 46 is also operable 
to track the satellite message data modulated on top of 
the Gold code ranging signal transmitted by each of the 
satellites 12 and to track the transition times of the 
satellite message data bit boundaries. 

10 The satellites 12 comprise satellites of a 

navigational satellite constellation that are in view of 
the ranging receiver 14 and/or other receivers 48 and 60 
at a particular time and that are operable to transmit 
positioning signals, such as navigation and ranging 

15 signals. It will be understood that the in-view 

satellites 12 of the satellite constellation continually 
change due to orbital movement of the satellites and 
rotational movement of the Earth. According to one 
embodiment, the system 10 comprises at least four 

20 satellites 12 at any particular time. 

The ranging receiver 14, the service center 18, and 
the geolocation processor 24 may comprise logic encoded 
in media. The logic comprises functional instructions 
for carrying out program tasks. The media comprises 

25 computer disks or other computer- readable media, 
application- specif ic integrated circuits, field- 
programmable gate arrays, digital signal processors, 
other suitable specific or general -purpose processors, 
transmission media or other suitable media in which logic 

30 may be encoded and utilized. 

In operation for a particular embodiment, the 
service center 18 provides data communication between the 
ranging receiver 14 and the geolocation processor 24 to 
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allow the ranging receiver 14 to generate a pseudorange 
estimate and the geolocation processor 24 to generate a 
geolocation. In addition, the service center 18 provides 
voice communication between a user of the wireless device 
5 16 and an operator at the service center 18 . The service 
center 18 may manage emergency services, vehicle - 
dispatching services, or other suitable services which 
utilize geolocation processing to locate wireless devices 
16. An operator at the service center 18 monitors and 

10 manages emergency calls, communications with vehicles via 
wireless devices 16, or other suitable communications. 

The receiver 48 receives positioning signals from 
the satellites 12 and processes these signals. The 
workstation 46 determines the location of a wireless 

15 device 16 based on the signals received at the receiver 
4 8 and based on information exchanged with the ranging 
receiver 14 via the wireless device 16 and the service 
center 18. 

The telephone handset 44 allows an operator at 
20 service center 18 to call the wireless device 16 or to 
respond to an incoming call from the wireless device 16. 
Calls to the wireless device 16 may also be initiated by 
the workstation 46 via a modem 40. Each modem 40 
provides data exchange capability with the ranging 
25 receiver 14 via the wireless device 16 over the same 
voice channel used by a caller to speak to the operator 
at the service center 18. According to one embodiment, 
each modem 4 0 also provides time-exchange facility with 
the ranging receiver 14 over the voice channel . 
30 The process may begin by an exchange of time signals 

between the ranging receiver 14 and the geolocation 
processor 24 to establish a real time at the ranging 
receiver 14 . This time is used to coordinate the use of 

DAL01:7560S9.1 



ATTORNEY'S DOCKET 
068796.0103 



20 



PATENT APPLICATION 



further time-critical information to be exchanged between 
the ranging receiver 14 and the geolocation processor 24 . 
While this time exchange is taking place between the 
ranging receiver 14 and the geolocation processor 24, the 
5 ranging receiver 14 may begin collecting and decimating 
satellite signal samples received from the satellites 12 
for a carrier frequency estimating process. It will be 
understood that, in accordance with one embodiment of the 
present invention, the ranging receiver 14 may be 
10 operable to determine real time for geolocation 
processing from the positioning signals without 
exchanging time signals with the geolocation px-ocessor 
24 . 

A direct method for quickly estimating the received 

15 carrier frequency for each of the satellites 12, 
including those of which the ranging receiver 14 has an 
obstructed view, may be used. According to one 

embodiment, this method uses signal squaring to remove 
the biphase modulation and uses narrow-band filtering to 

20 reduce the noise. In one embodiment, at least two bits 
per signal sample and a final stage of down- conversion to 
near baseband by a heterodyne step are used. 

The received carrier frequency estimates include the 
down -conversion frequency uncertainty due to local 

25 oscillator manufacturing tolerances, thermal drift and/or 
other conditions, thereby providing a substantially exact 
frequency that may be applied to compensate for the 
Doppler frequency shift and these other errors, even over 
the long periods of signal averaging needed to reach the 

30 required sensitivity for successful signal detection. 

In one embodiment of the invention, this direct 
determination of the received carrier frequencies to high 
resolution (within about one Hz for a one second sampling 
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of the signals received from the satellites 12) allows 
the identity of satellites 12 contributing the specific 
pattern of Doppler- frequency differences, which are equal 
to the carrier frequency offset (CFO) differences, to be 
5 determined by a pattern match between the Doppler 
differences computed for the current positions of the 
satellites in the constellation and target matching 
patches on the surface of the Earth. This may be 
accomplished by sending the CFO results to the 

10 geblocation processor 24 as soon as they are determined 
by the ranging receiver 14, along with the time 
corresponding to the beginning of the collection of the 
signal samples from which the CFOs were determined. The 
ranging receiver 14 then continues collecting and storing 

15 samples of the signal that are suitable for estimating 
pseudorange . 

While the signal samples for pseudorange estimation 
are being collected, the geolocation processor 24 
identifies the satellites 12 that are in view at the 

20 unknown location and, for each of these satellites 12, 
extracts that portion of the satellite message data that 
is being transmitted during sample collection. The 
geolocation processor 24 may also identify the actual 
Doppler frequency shift frequency to be used in the time- 

25 dilation correction, as described in FIGURE 13C. Time- 
dilation is proportional to the total accumulated Doppler 
frequency shift. The geolocation processor 24 then 
transmits to the ranging receiver 14 a satellite ID for 
each identified satellite 12 to enable rapid computation 

30 of pseudorange without searching through Doppler 
frequencies and satellite codes and transmits the 
satellite message data fragment to enable removal of the 
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data modulation from the signals so that linear summation 
may be performed for increased processing gain. 

Also while the ranging receiver 14 is collecting 
signal samples for pseudorange estimation, the ranging 
5 receiver 14 begins processing the signals for the 
pseudorange content . 

FIGURE 2 is a block diagram illustrating the ranging 
receiver 14 in accordance with one embodiment of the 
present invention. The ranging receiver 14 is a mobile 

10 device that may be easily transported from one location 
to another. In the illustrated embodiment, the ranging 
receiver 14 is coupled to a wireless device 16 that 
comprises a digital cellular telephone. However, as 
described in more detail above in connection with 

15 FIGURE 1, the wireless device 16 may comprise any 
suitable device capable of wireless communication. 

The ranging receiver 14 comprises an antenna 2 02 
that is capable of receiving positioning signals. 
According to one embodiment, the antenna 202 comprises a 

20 GPS antenna that is capable of receiving GPS signals. 
For this embodiment, the antenna 202 may comprise a 
circularly polarized ceramic patch antenna or other 
suitable antenna. The antenna 202 is coupled to a filter 
and low-noise amplifier 204. 

25 In one embodiment, the filter and low-noise 

amplifier 204 is operable to limit bandwidth to 
approximately 2 MHz (centered at 1575.42 MHz). The 
filter and low-noise amplifier 204 is coupled to a down- 
converter 206. The down-converter 206 is operable to 

3 0 down- convert the signal to baseband. According to one 
embodiment, the down -converter is operable to down- 
convert the signal using multiple stages of 
amplification, filtering and super- heterodyning. The 
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down-converter 206 may comprise a radio frequency/ 
intermediate frequency (RF/IF) down -converter or other 
suitable down- converter . 

An automatic gain control 208 is coupled to the 
down- converter 206 and to an amplifier 210. The 
automatic gain control 208 controls the stages of 
amplification for the down- converter 206. The amplifier 
210, which is also coupled to the down- converter 206, 
amplifies the signal from the down -converter 206 based on 
a control signal from the automatic gain control 208. 

The amplifier 210 may be operable to limit bandwidth 
to approximately 1 MHz or to another suitable bandwidth 
less than 2 MHz . In this way, the signal-to-noise ratio 
may be improved as compared to a signal with a bandwidth 
of approximately 2 MHz. For example, for the embodiment 
in which the bandwidth is limited to approximately 1 MHz, 
a gain of about 3dB made be obtained with a loss of about 
50% of the noise and only about 10% of the signal. 

An analog-to-digital converter 212 is coupled to the 
amplifier 210 and to the automatic gain control 208. The 
analog-to-digital converter 212 comprises a symmetrical 
multi-bit quantizer that is operable to digitize analog 
signals from the down- converter 206. The analog-to- 
digital converter 212 is symmetrical in order to minimize 
the generation of unnecessary harmonic products that may 
cause alias products in the band of the desired signal. 

For example, the analog- to-digital converter 212 may 
digitize the signals at a minimum of 2.048 
Megasamples/second and at a quantizer resolution of a 
minimum of two bits. Alternatively, the analog-to- 
digital converter 212 may digitize the signals at a 
minimum of 4.096 or 8.192 Megasamples/second, in order to 
improve pseudorange timing resolution, and then further 
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into four or more bits, in-phase and quadrature -phase 
(I/Q) samples, in order to improve both frequency and 
amplitude resolution and to minimize the generation of 
harmonic products due to the inherent non-linearity of 
5 the quantizing process. 

Symmetry may be ensured by biasing the quantizer 
levels such that zero crossing takes place in the center 
of the first quantizer level, allowing equal statistical 
representation of the signal in the positive and negative 

10 segments. This is especially useful for small numbers of 
bits of resolution. 

An output of the analog-to-digital converter 212 
controls the automatic gain control 208, which feeds back 
to the amplifier 210 and the down-converter 206 as 

15 described above. This may be used to ensure that, on 
average, the quantized signal reaches maximum quartile 
less than 30% of the time, which is consistent with the 
expected statistics of noise-like signals. However, it 
will be understood that this output of the analog-to- 

20 digital converter 212 may be used to ensure that the 
quantized signal reaches maximum quartile less than any 
suitable percentage of the time. 

The ranging receiver 14 also comprises a local 
oscillator and timing generator 214 that is coupled to 

25 the down-converter 206 and to the analog-to-digital 
converter 212 . The local oscillator and timing generator 
214 may comprise a temperature -compensated crystal 
oscillator (TCXO) , a digitally-compensated crystal 
oscillator (DCXO) , or other suitable oscillator operable 

30 to drive a digital device. Various timing signals may be 
derived from the local oscillator and timing generator 
214, such as frequencies for the various stages of down- 
conversion for the down -converter 206, the sampling clock 
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for the analog- to-digital converter 212, and other 
suitable timing signals. 

A signal sample memory 216 is coupled to the analog- 
to-digital converter 212. The signal sample memory 216 
5 is operable to temporarily store digital signal samples 
from the analog-to-digital converter 212. For example, 
the signal sample memory 216 may comprise a random access 
memory (RAM) or other suitable data store. 

The ranging receiver 14 also comprises a digital 

10 signal processor 218 that is coupled to the signal sample 
memory 216. The digital signal processor 218 comprises a 
processor operable to process stored instructions. For 
example, the digital signal processor 218 may comprise a 
Motorola DSP 56654 manufactured by MOTOROLA, INC. of 

15 Schaumburg, Illinois, a Texas Instruments TMS 320VCSSIO 
manufactured by TEXAS INSTRUMENTS of Dallas, Texas, an 
ADSP 21161 SHARC DSP manufactured by ANALOG DEVICES of 
Norwood, Massachusetts, or other suitable processor. 

A receiver memory 220 is coupled to the digital 

20 signal processor 218. The receiver memory 220 is 
operable to store instructions and other data for the 
digital signal processor 218. The receiver memory 220 
may comprise a plurality of program and working memories, 
which may comprise any suitable combination of read-only 

25 memories (ROMs), RAMs, and other suitable data stores. 
It will be understood that the receiver memory 220 may 
also comprise the signal sample memory 216 . 

A modem function 222 is coupled to the digital 
signal processor 218 and is operable to provide an 

3 0 interface through which communication may be provided 
between the digital signal processor 218 and the wireless 
device 16 and ultimately with a geolocation processor 24 
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through the wireless carrier network 20, the POTS network 
22, and a service center 18. 

Signal lines 224 and 226 couple the wireless device 
, 16 to the modem function 222 and to the digital signal 
5 processor 218. According to the illustrated embodiment, 
the signal lines 224 and 226 are coupled to an external 
component 228 of the wireless device. It will be 
understood, however, that the signal lines 224 and 226 
may be otherwise coupled to the wireless device 16 based 

10 on which embodiment of the ranging receiver 14 is 
implemented, as described above in connection with 
FIGURE 1. For example, for the embodiment in which the 
ranging receiver 14 is an integral component of the 
wireless device 16, the signal lines 224 and 226 may 

15 comprise internal connections within the wireless device 
16 which need not be coupled to an external component 
228 . 

The ranging receiver 14 also comprises power 
management circuitry 23 0 that is coupled to the digital 

20 signal processor 218. The power management circuitry 230 
may comprise one or more batteries or other suitable 
power sources. The power management circuitry 230 is 
operable to provide power to the ranging receiver 14 and 
to control the availability of power on the basis of 

25 control signals on the signal lines 224 and 226 from the 
wireless device 16 and/or other suitable control signals. 
The ranging receiver 14 also comprises a real-time clock 
232 that is operable to be calibrated to a remote clock 
across a wireless link or based on satellite message 

30 data. 

FIGURE 3 is a block diagram illustrating the ranging 
receiver 14 in accordance with another embodiment of the 
present invention. In the illustrated embodiment, an 
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integrated device 250 comprises the function of a ranging 
receiver 14 integrated with a wireless device 16 that 
comprises a digital cellular telephone. However, as 
described in more detail above in connection with 
5 FIGURE 1, the wireless device 16 may comprise any 
suitable device capable of wireless communication. 

The integrated device 250 comprises an antenna 252 
that is capable of receiving positioning signals. 
According to one embodiment, the antenna 252 comprises a 

10 GPS antenna that is capable of receiving GPS signals from 
the satellites 12. For this embodiment, the antenna 252 
may comprise a circularly polarized ceramic patch antenna 
or other suitable antenna. The integrated device 250 
also comprises an antenna 254 that is operable to 

15 communicate with a base station 25S. Thus, for the 
digital cellular telephone embodiment, the antenna 254 is 
capable of receiving and transmitting cellular 
communication signals. 

The antenna 2 52 is coupled to a positioning signal 

20 (^^PS") preselect filter 260, and the antenna 254 ■ is 
coupled to a cellular signal preselect filter 262. The 
output from each of the filters 260 and 262 is coupled to 
a signal selection switch 264. 

The signal selection switch 264 is operable to 

25 select either the positioning signals or the cellular 
communication signals for processing by the integrated 
device 250. The output of the signal selection switch 
264, which corresponds to the output of the filter 260 or 
262 selected by the switch 264, is coupled to additional 

30 cell phone circuitry 266. The additional cell phone 
circuitry 266 is operable to process positioning signals 
as described above in connection with FIGURE 2 and to 
process cellular communication signals in a manner 
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similar to a standard wireless device 16 without an 
integrated ranging receiver 14. 

FIGURE 4 is a flow diagram illustrating a method for 
obtaining an initial calibration location estimate based 
on an initial trigger for the ranging receiver 14 in 
accordance with one embodiment of the present invention. 
The method begins at step 400 where the ranging receiver 
14 detects an initial trigger. 

A trigger may comprise a request generated by the 
wireless device 16 automatically, a request generated by 
a user of the wireless device 16 such as by pushing a 
button 3 0 on the ranging receiver 14, a request from an 
external source such as the Internet, a geolocation 
processor 24 or other suitable external source, or any 
other suitable trigger. The initial trigger detected in 
step 400 may comprise, for example, a trigger that is 
automatically generated in response to turning on the 
wireless device 16 that is coupled to the ranging 
receiver 14. In addition, as described in more detail 
below, the initial trigger may comprise a trigger that is 
automatically generated after a previous initial trigger 
was followed by an unsuccessful attempt to obtain a 
calibration location estimate. 

At decisional step 402, a determination is made 
regarding the status of a superframe in the ranging 
receiver 14. According to one embodiment, the superframe 
comprises 25 frames, each of which further comprises five 
subframes. The superframe comprises 37,500 bits, which 
are transmitted from the satellites 12 at 50 bits per 
second. Each of the 25 frames comprises a 1500-bit 
message that includes timing, identification, and 
ephemeris data, along with portions of the 



DALOl: 756059.1 



ATTORNEY'S DOCKET 
068796 .0103 



29 



PATENT APPLICATION 



constellation's almanac data and various telemetry, 
checksum, and counterwords for each satellite 12. 

Changes in the message relating to updates to the 
ephemeris and almanac data occur at infrequent intervals, 
usually between four and 24 hours, and affect only a 
small portion of the message for each update. Therefore, 
if the superframe in the ranging receiver 14 is 
relatively current, the message in that superframe may be 
used by the ranging receiver 14. 

If no superframe is present in the ranging receiver 
14 or if the superframe is expired, the method follows 
the Expired branch from decisional step 402 to step 404. 
The superframe may be considered expired if the 
superframe is greater than 96 hours old. However, it 
will be understood that the superframe may be considered 
expired based on another suitable amount of time or other 
suitable criteria without departing from the scope of the 
present invention. At step 404, the ranging receiver 14 
downloads the superframe from the satellites 12 or from 
the geolocation processor 24 . 

If the superframe in the ranging receiver 14 is not 
current, the method follows the Not Current branch from 
decisional step 402 to step 406. The superframe may be 
considered not current if the superframe is between 4 and 
96 hours old. However, it will be understood that the 
superframe may be considered not current based on another 
suitable amount of time or other suitable criteria 
without departing from the scope of the present 
invention. 

At Step 406, the ranging receiver 14 downloads an 
update to the superframe from the satellites 12 or from 
the geolocation processor 24 . Any data obtained from the 
geolocation processor 24 in steps 404 or 406 may be 
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obtained over the voice channel using the modems 4 0 or 
via an alternate pathway using, for example, a control 
channel for the wireless carrier network 20 or other 
suitable data pathway. 
5 Returning to decisional step 402, if the superframe 

in the ranging receiver 14 is current, the method follows 
the Current branch from decisional step 402 to decisional 
step 408. Also, after a superframe has been downloaded 
in step 404 or an update has been downloaded in step 406, 

10 the method continues to decisional step 408. The 
superframe may be considered current if the superframe is 
less than 4 hours old. However, it will be understood 
that the superframe may be considered current based on 
another suitable amount of time or other suitable 

15 criteria without departing from the scope of the present 
invention. 

At decisional step 408, a determination is made 
regarding whether or not a fresh calibration location 
estimate is available. According to one embodiment, a 

20 fresh calibration location estimate comprises a location 
estimate that has been performed within the past ten 
minutes. However, it will be understood that a 

calibration location estimate may be considered fresh 
after any suitable period of time or based on other 

25 suitable criteria without departing from the scope of the 
present invention. 

If no fresh calibration location estimate is 
available, the method follows the No branch from 
decisional step 408 to step 410. At step 410, the 

3 0 location circuitry of the ranging receiver 14 is turned 
on. Thus, for example, the power management circuitry 
230 may provide full power to each of the components 204, 
206, 208, 210, 212, 214, 216, 218, 220 and 222. At step 
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412, the ranging receiver 14 begins collecting samples 
for' CFO extraction and attempts to perform CFO extraction 
based on the collected samples. 

At decisional step 414, a determination is made 
regarding whether or not the CFO extraction was 
successful. If the CFO extraction was successful, the 
method follows the Yes branch from decisional step 414 to 
step 416. At step 416, the ranging receiver 14 enters a 
stand-alone mode for pseudorange estimation. At step 
418, the ranging receiver 14 estimates pseudorange. 
According to one embodiment, the ranging receiver 14 
estimates pseudorange in accordance with the method of 
FIGURES 7 and 8 . 

At decisional step 420, a determination is made 
regarding whether or not the real-time clock 232 has the 
current time. If the real-time clock 232 does not have 
the current time, the method follows the No branch from 
decisional step 420 to step 422. 

At step 422, the real-time clock 232 for the ranging 
receiver 14 is calibrated, or synchronized, to a 
standardized time of day (TOD) , such as GPS time, across 
the wireless link or to the TOD information received in 
the positioning signals. Once calibrated, the real-time 
clock 232 keeps track of real time so that the timing of 
events in the ranging receiver 14, such as message data 
bit transition times, can be accurately tracked. From 
step 422, the method continues to step 424. 

Returning to decisional step 420, if the real-time 
clock 232 has the current time, the method follows the 
Yes branch from decisional step 420 to step 424. At step 
424, the ranging receiver 14 performs a calibration 
location estimate. 
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At step 426, the ranging receiver 14 determines its 
receiver bias based on geometric and measured data. The 
geometric data may comprise satellite message data, 
approximate location data and/or other suitable data 
5 useful for estimating Doppler frequency shift. The 
approximate location data may comprise a previous 
location estimate, such as a calibration location 
estimate, or general location information provided by an 
external source, such as a cell site identification 

10 provided by a PSAP or an estimate provided by a user of 
the wireless device 16. According to one embodiment, the 
ranging receiver 14 determines its receiver bias by first 
determining the actual Doppler frequency shift for the 
satellite signals being processed based on the geometric 

15 data from each satellite 12. Based on the actual Doppler 
frequency shifts (estimated from the geometric data) and 
the measured CFOs, the ranging receiver 14 calculates the 
receiver frequency bias for the ranging receiver 14, 
which is stored for future use. At step 428, the ranging 

20 receiver 14 reverts to the low-power standby condition 
and waits for another trigger. 

At step 430, the ranging receiver 14 sets a timer 
for a recalibration trigger. Thus, after the timer 
expires, a recalibration trigger will be generated in 

25 step 434 to prompt the ranging receiver 14 to perform 
another location estimate. The method then returns to 
step 43 0 where the timer is reset for another 
recalibration trigger. In this way, the location 

estimate may be automatically updated periodically. 

3 0 According to one embodiment, the timer set in step 

430 expires after approximately 10 minutes. However, it 
will be understood that the timer may expire after any 
suitable amount of time without departing from the scope 
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of the present invention. For example, according to one 
embodiment, the timer may expire after a varying amount 
of time based on comparing consecutive location 
estimates. Thus, for this embodiment, the timer may 
expire after a reduced amount of time when consecutive 
location estimates indicate that the ranging receiver 14 
is moving at greater than a certain velocity and may 
expire after an increased amount of time when consecutive 
location estimates indicate that the ranging receiver 14 
is moving at less than a certain velocity. 

Returning to decisional step 414, if the CFO 
extraction was not successful, the method follows the No 
branch from decisional step 414 to step 440. At step 
44 0, the ranging receiver 14 sets a timer for an initial 
trigger. Thus, after a specified amount of time has 
passed (i.e., after the timer expires), an initial 
trigger will be generated in step 442. The method then 
returns to step 400 where the ranging receiver 14 detects 
the initial trigger and the method begins again. In this 
way, the calibration location estimate may be re- 
attlmpted automatically at a later time. According to 
one embodiment, the specified amount of time is 
approximately 10 minutes. However, it will be understood 
that the specified amount of time may be any suitable 
i length of time without departing from the scope of the 
present invention. 

Returning to decisional step 408, if a fresh 
calibration location estimate is determined to be 
available, the method follows the Yes branch from 
D decisional step 408 to step 430 where the timer is set 
for a recalibration trigger as described above. In this 
way, the calibration location estimate may be updated 
automatically after the specified period of time. 
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FIGURE 5 is a flow diagram illustrating a method for 
obtaining a location estimate based on a non- initial 
trigger for the ranging receiver 14 in accordance with 
one embodiment of the present invention. The method 
begins at step 500 where the ranging receiver 14 detects 
a non- initial trigger. 

According to one embodiment, the non-initial trigger 
may comprise a recalibration trigger that is 
automatically generated by the ranging receiver 14 
periodically based on a timer. Thus, when the 

recalibration trigger is generated, the timer for the 
recalibration trigger may be automatically reset to 
trigger another location estimate after a specified 
amount of time has passed. The non- initial trigger may 
also comprise a trigger that is generated by a user of 
the wireless device 16 pushing a button 3 0 on the ranging 
receiver 14, generated by a service requester 62, or 
based on any other suitable request. 

At decisional step 502, a determination is made 
regarding the status of a superframe in the ranging 
receiver 14. If no superframe is present in the ranging 
receiver 14 or if the superframe is expired, the method 
follows the Expired branch from decisional step 502 to 
step 504. The superframe may be considered expired if 
the superframe is greater than 96 hours old. However, it 
will be understood that the superframe may be considered 
expired based on another suitable amount of time or other 
suitable criteria without departing from the scope of the 
present invention. At step 504, the ranging receiver 14 
downloads the superframe from the satellites 12 or from 
the geolocation processor 24. 

If the superframe in the ranging receiver 14 is not 
current, the method follows the Not Current branch from 
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decisional step 502 to step 506. The superframe may be 
considered not current if the superframe is between 4 and 
96 hours old. However, it will be understood that the 
superframe may be considered not current based on another 
5 suitable amount of time or other suitable criteria 
without departing from the scope of the present 
invention. At step 506, the ranging receiver 14 

downloads an update to the superframe from the satellites 
12 or from the geolocation processor 24. 

10 Any data obtained from the geolocation processor 24 

in steps 504 or 506 may be obtained over the voice 
channel using the modems 4 0 or via an alternate pathway 
using, for example, a control channel for the wireless 
carrier network 2 0 or other suitable data pathway. 

15 Returning to decisional step 502, if the superframe 

in the ranging receiver 14 is current, the method follows 
the Current branch from decisional step 502 to step 508. 
Also, after a superframe has been downloaded in step 504 
or an update has been downloaded in step 506, the method 

20 continues to step 508. The superframe may be considered 
current if the superframe is less than 4 hours old. 
However, it will be understood that the superframe may be 
considered current based on another suitable amount of 
time or other suitable criteria without departing from 

25 the scope of the present invention. 

At step 508, the location circuitry of the ranging 
receiver 14 is turned on. At step 510, the ranging 
receiver 14 begins collecting samples for CFG extraction 
and attempts to perform CFG extraction based on the 

30 collected samples. 

At decisional step 512, a determination is made 
regarding whether or not the CFG extraction was 
successful. If the CFO extraction was successful, the 

DALOl: 756059.1 



ATTORNEY'S DOCKET 
068796 .0103 



36 



PATENT APPLICATION 



method follows the Yes branch from decisional step 512 to 
step 514. At step 514, the ranging receiver enters the 
stand-alone mode for pseudorange estimation. At step 
516, the ranging receiver 14 estimates pseudorange. 
5 According to one embodiment, the ranging receiver 14 
estimates pseudorange in accordance with the method of 
FIGURES 7 and 8. 

At decisional step 518, a determination is made 
regarding whether or not the real-time clock 232 has the 

10 current time. If the real-time clock 232 does not have 
the current time, the method follows the No branch from 
decisional step 518 to step 520. 

At step 520, the real-time clock 232 for the ranging 
receiver 14 is calibrated, or synchronized, to a remote 

15 clock across the wireless link or based on TOD 
information in the positioning signals. Once calibrated, 
the real-time clock 232 keeps track of real time so that 
the timing of events in the ranging receiver 14, such as 
message bit transition times, can be accurately tracked. 

20 From step 520, the method continues to step 522. 

Returning to decisional step 518, if the real-time 
clock 232 has the current time, the method follows the 
Yes branch from decisional step 518 to step 522. At step 
522, the ranging receiver 14 performs the location 

25 estimate. 

At step 524, the ranging receiver 14 determines its 
receiver bias based on geometric data. The geometric 
data may comprise satellite message data, approximate 
location data and/or other suitable data useful for 
3 0 estimating CFO for compensation during pseudorange 
processing in the geometric mode. The approximate 
location data may comprise a previous location estimate, 
such as a calibration location estimate, or general 
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location information provided by an external source, such 
as a cell site identification provided by a PSAP or an 
estimate provided by a user of the wireless device 16. 
According to one embodiment, the ranging receiver 14 
5 determines its receiver bias by first calculating the 
actual Doppler frequency shift for each of the satellite 
signals being received based on the geometric data for 
each satellite 12 and the most recent location estimate. 
Based on these actual Doppler frequency shifts, the 

10 ranging receiver 14 calculates the difference between the 
actual measured CFOs and the calculated Doppler frequency 
shifts to determine the frequency bias for the ranging 
receiver 14, which is stored for future use. From step 
524, the method continues to step 526. 

15 Returning to decisional step 512, if the CFG 

extraction was not successful, the method follows the No 
branch from decisional step 512 to step 528 . At step 
528, the ranging receiver 14 enters the geometric mode. 
Thus, when the ranging receiver 14 can successfully 

20 perform CFO extraction, the ranging receiver 14 enters 
the stand-alone mode in order to determine a pseudorange. 
However, when the ranging receiver 14 is unable to 
perform CFO extraction due to the weak level of the 
signals being received at the ranging receiver 14, the 

25 ranging receiver 14 enters the geometric mode in order to 
determine a pseudorange . According to one embodiment , 
the ranging receiver 14 enters the geometric mode if 
attempted CFO extraction is unsuccessful within a pre- 
defined period of time. Depending on the application or 

30 embodiment, the pre-defined period of time may comprise 
approximately 1 to 3 0 seconds or any other suitable 
period of time. 
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At step 53 0, the ranging receiver 14 obtains 
geometric data for the location estimating process. As 
used herein "to obtain data" means to retrieve or receive 
data from a local memory store or from a remote device. 
5 According to one embodiment, the geometric data comprises 
satellite message data and approximate location data. 
The approximate location data may comprise a previous 
location estimate, such as a calibration location 
estimate, or general location information provided by an 

10 external source, such as a cell site identification 
provided by a PSAP or an estimate provided by a user of 
the wireless device 16 . 

At step 532, the ranging receiver 14 estimates the 
CFO for each satellite 12. According to one embodiment, 

15 these CFOs are estimated by summing the Doppler frequency 
shifts (which may be determined based on the geometric 
data) and the receiver frequency bias (which was 
previously determined during the calibration process) . 

At step 534, the ranging receiver 14 obtains message 

20 data fragments from the geolocation processor 24 or from 
the results of the method for predicting satellite 
message data described in FIGURE 12 or other suitable 
prediction method. At step 536, the ranging receiver 14 
estimates pseudorange . According to one embodiment, the 

25 ranging receiver 14 estimates pseudorange in accordance 
with the method of FIGURES 13A, 13B and 13C. At step 
538, the ranging receiver 14 performs the location 
estimate . 

From steps 524 and 538, the method continues to step 
30 526 where the ranging receiver 14 stores the location 
estimate, receiver bias and current time. At step 540, 
the ranging receiver 14 reverts to the low-power standby 
condition. At decisional step 542, a determination is 

DALOl: 756059.1 



ATTORNEY'S DOCKET 
068796.0103 



39 



PATENT APPLICATION 



made regarding whether the trigger detected in step 500 
was a recalibration trigger. If the trigger was not a 
recalibration trigger, the method follows the No branch 
from decisional step 542 to step 544. At step 544, the 
5 ranging receiver 14 sends the location estimate to the 
requester who generated the trigger, at which point the 
method comes to an end. However, if the trigger was 
determined to be a recalibration trigger in decisional 
step 542, the method follows the Yes branch from 

10 decisional step 542 and comes to an end. 

FIGURES 6A and 6B are flow diagrams illustrating a 
method for direct CFO extraction in the system 10 in 
accordance with one embodiment of the present invention. 
The method begins at step 600 where the location 

15 circuitry of the ranging receiver 14 is turned on. At 
step 6 02, a decimated output matrix, DM, of dimension 
(P,W) with element pointer p,w, is initialized and the 
decimated output pointer, p, is initialized to 1. 

At step 604, sample sizes SI and S2 are initialized, 

20 a maximum number of samples, MAX, is initialized to SI, 
and a sample number, s, is set to 1. Also at step 604, a 
rate correction indicator, RC, is set to 0 . In one 
embodiment, the sample size SI is set to correspond to a 
sampled duration of approximately one second and the 

25 sample size S2 is set to correspond to a sampled duration 
between approximately 5 and 100 seconds. However, it 
will be understood that the sample sizes may be set to 
any suitable number without departing from the scope of 
the present invention. 

30 At step 6 06, a number W of rate values, R„, is 

initialized and w is set to a value of one. In one 
embodiment, the number of values is selected to ensure 
that the error range between values is less than 
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0.1 Hz/second. The number of values may be selected as 
the maximum rate range divided by the error range. 
According to one embodiment, the maximum rate range may 
reach approximately 1.0 Hz/second and the number of 
5 values may be selected to be 1/(0.1 - -0.1), or 5. It 
will be understood, however, that any suitable criteria 
may be used in selecting the number of values. 

At step 608, a decimation counter is initialized. 
According to one embodiment, the decimation counter is 

10 initialized to a decimation rate, DR, that corresponds to 
the rate at which the samples will be decimated. This 
decimation rate may be chosen in relation to the signal - 
sampling rate {which may be 1,024 kHz x 2"^, with n = 1, 2 
or 3 , or other suitable rate) . For example, the signal 

15 may be sampled directly to result in down- conversion to 
baseband at these sampling rates or at higher rates to 
produce results that may be easily converted to I/Q 
samples. According to one embodiment, the decimation 
rate may be chosen to yield a signal with at least a 20 

20 kHz bandwidth, thereby allowing for the maximum range of 
the CFOs and the doubling in frequency caused by the 
squaring step. Thus, for a 2,048 kHz complex sampling 
rate (n=l, above), the decimation rate (DR) may be chosen 
to be approximately 100. 

25 At step 610, the ranging receiver 14 collects a new 

signal sample. At step 611, the ranging receiver 14 
increments the sample counter, s. At step 612, the 
ranging receiver 14 squares the signal sample in order to 
recover the original carrier on which the biphase Gold 

30 code is modulated and to double its original frequency. 
At step 614, the ranging receiver 14 adds the squared 
sample value to an accumulator. According to one 
embodiment, the squared sample is multiplied by a 
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filtering coefficient identified by the decimation 
counter before being added to the accumulator. At step 
616, the ranging receiver 14 decrements the decimation 
counter . 

5 At decisional step 618, a determination is made 

regarding whether or not the decimation counter is equal 
to zero. If the decimation counter is not equal to zero, 
the method follows the No branch from decisional step 618 
and returns to step 610 for the collection of another 

10 signal sample. However, if the decimation counter is 
equal to zero, the method follows the Yes branch from 
decisional step 618 to decisional step 620. 

At decisional step 620, a determination is made 
regarding whether or not RC equals 1 . According to the 

15 illustrated embodiment, rate correction is applied when 
RC equals 1 and is not applied when RC equals 0 . 
However, it will be understood that rate correction may 
be applied based on any suitable vales for RC or based on 
other suitable criteria without departing from the scope 

20 of the present invention. For the illustrated 

embodiment, if RC equals 1, the method follows the Yes 
branch from decisional step 620 to step 621. At step 
621, w is set to a value of one. 

At step 622, the ranging receiver 14 gets the w*^*^ 

25 Doppler correction. In order to determine the w*^*^ Doppler 
correction, the ranging receiver 14 computes a pointer to 
a table of sine and cosine values. These tabular values 
are used to implement multiplication of the samples by 
the phase rotation vector e"^'"*^, i.e., cos (cot ) - j • sin (cot ) , 

3 0 where j represents the square root of -1 and co is the 
frequency term 27t{D„+Rwt), where D„ is the current, or w*"*", 
Doppler frequency shift and R„ is the current, or w'^*', 
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Doppler rate. Since D„ is unknown, it is set to zero, and 
since the estimate producing the largest response will 
identify the appropriate rate Rw, D„ will come out of the 
spectral analysis of each of the W rate -compensated 
5 arrays as the CFO. The variable t is the time from the 
beginning of the analysis to the center of the current 
decimated sample, i.e., t = i/f sampie-rate/ and i = signal 
sample number at the center of this decimated sample. 

At step 624, the ranging receiver 14 complex 

10 multiplies the complex decimated sample (I/Q sampling) 
accumulated at step 614 by the w*^"^ complex Doppler 
correction. At step 626, the ranging receiver 14 stores 
the results of the complex multiplication at a location 
identified by p,w in the decimated output matrix in the 

15 signal sample memory 216. At step 628, the ranging 
receiver 14 increments w. 

At decisional step 63 0, a determination is made 
regarding whether or not w is greater than W. If w is 
not greater than W, the method follows the No branch from 

20 decisional step 630 and returns to step 622, where the 
ranging receiver 14 gets the Doppler correction for the 
incremented w. However, if w is greater than W, the 
method follows the Yes branch from decisional step 63 0 to 
step 632. At step 632, the ranging receiver 14 

25 increments p. 

Returning to decisional step 620, if RC does not 
equal 1, the method follows the No branch from decisional 
step 620 to step 634. At step 634, the value in the 
accumulator described in step. 614 is stored at a location 

30 identified by p,w in the signal sample memory 216. At 
this point, the method continues to step 632 where the 
ranging receiver 14 increments p. From step 632, the 
method continues to decisional step 635. 
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At decisional step 635, a determination is made 
regarding whether or not p is equal to the maximum number 
of decimated samples, MAX/DR. If P is not equal to 
MAX/DR, the method follows the No branch from decisional 
step 635 and returns to step 608 where the decimation 
counter is re- initialized for the process of collecting 
additional signal samples. However, if P is equal to 
MAX/DR, the method follows the Yes branch from decisional 
step 635 to step 636. 

At step 636, the ranging receiver 14 attempts to 
find significant carriers in the sample sets. According 
to one embodiment, the ranging receiver 14 attempts to 
find significant carriers in accordance with the method 
of FIGURE 6B. While the ranging receiver is attempting 
to find significant carriers, the method continues to 
decisional step 640. 

At decisional step 640, a determination is made 
regarding whether or not s is greater than MAX. If s is 
not greater than MAX, the method follows the No branch 
from decisional step 640 and returns to step 608 where 
the decimation counter is re-initialized. However, if s 
is greater than MAX, the method follows the Yes branch 
from decisional step 640 to step 642. At step 642, RC is 
set to a value of one, and MAX is increased from SI to 
i the greater value S2 . 

At decisional step 644, a determination is made 
regarding whether or not s is greater than MAX. If s is 
not greater than MAX, the method follows the No branch 
from decisional step 644 and returns to step 608 where 
3 the decimation counter is re-initialized. However, if s 
is greater than S2 , the method follows the Yes branch 
from decisional step 644 to step 646. At step 646, a 
flag is set to indicate that the signal is too weak for 
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successful CFO extraction, at which point the method 
comes to an end. 

FIGURE 6B is a flow diagram illustrating a method 
for attempting to find significant carriers in the method 
5 of FIGURE 6A in accordance with one embodiment of the 
present invention. The method begins at step 650 where 
the decimated sample sets stored in the signal sample 
memory 216, each of which corresponds to a Doppler-rate 
compensation estimate, are Fourier-transformed into the 

10 frequency domain. 

At step 652, the spectrum of each set is searched 
for significant carriers, i.e., carriers that are 
sufficiently above a noise-background threshold. 
According to one embodiment, a carrier is sufficiently 

15 above the noise-background threshold when the carrier is 
at least 6dB above the noise -background threshold. 
However, it will be understood that a carrier may be 
sufficiently above the noise -background threshold when 
the carrier is any suitable level above the noise- 

20 background threshold. 

In a particular embodiment, the six highest peaks 
are identified as possible carriers and the 20 next 
highest peaks are averaged to establish a noise- 
background threshold. The possible carriers may be 

25 confirmed as carriers that are sufficiently above the 
noise -background threshold when the carriers are at least 
four times, or the equivalent of 6dB, higher than the 
threshold. It will be understood that the noise- 
background threshold may be otherwise suitably determined 

3 0 without departing from the scope of the present 
invention. In addition, the carriers may be considered 
to be sufficiently above the noise -background threshold 
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when the carriers are at any suitable level above the 
threshold. 

At decisional step 654, a determination is made 
regarding whether a specified minimum number of 
5 significant carriers sufficiently above the noise- 
background threshold has been found. According to one 
embodiment, the specified minimum number is four. 
According to another embodiment, the specified minimum 
number may be greater than four since the availability of 

10 more pseudorange estimates improves the quality or 
confidence level of a location estimate. However, it 
will be understood that any suitable number of 
significant carriers may be specified. 

If the specified number of significant carriers has 

15 been found, the method follows the Yes branch from 
decisional step 654 to step 656. At step 656, the 
ranging receiver 14 stores the CFG results. At step 658, 
the direct CFO extraction of the method of FIGURE 6A is 
ended due to the successful attempt to find significant 

20 carriers. In this way, CFOs can be extracted from 
approximately one second of signal sample processing from 
signals obstructed up to approximately lOdB of 
attenuation or can be extracted from approximately 10 
seconds of signal sample processing from signals 

25 obstructed up to approximately 15dB of attenuation. 

Returning to decisional step 654, if the specified 
number of significant carriers has not been found, the 
method follows the No branch from decisional step 654 and 
returns to the method of FIGURE 6A for the continuation 

3 0 of direct CFO extraction. 

FIGURE 7 is a flow diagram illustrating a method for 
determining a pseudorange in the ranging receiver 14 in 
the stand-alone mode in accordance with one embodiment of 
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the present invention. The method begins at step 700 
where the ranging receiver 14 begins to collect samples 
for pseudorange estimation from the signals received from 
the satellites 12. At step 702, the ranging receiver 14 
5 selects one of the CFOs determined from the direct CFO 
extraction process of FIGURE 6 for the satellites 12. At 
step 704, the ranging receiver 14 begins to process the 
pseudorange samples by compensating the samples for the 
selected CFO. As used herein, ^'^pseudorange samples" 
10 means signal samples from which pseudorange may be 
estimated. 

At step 706, the ranging receiver 14 removes 
satellite message data modulation from the pseudorange 
samples by multiplying out the biphase message data 

15 modulation, as described in more detail below in 
connection with FIGURES 8A and SB. 

At step 708, the ranging receiver 14 linearly 
combines, or stacks, the samples into a single epoch 
duration result for each satellite 12. At step 710, the 

20 ranging receiver 14 compensates for time dilation. At 
step 712 , the ranging, receiver 14 correlates the Gold 
codes for the satellites 12 to give the pseudoranges from 
those satellites 12. According to one embodiment, the 
ranging receiver 14 correlates the Gold codes by Fourier 

25 transform means. For this embodiment, each CFO- 

compensated stack of samples which has been Fourier- 
transformed is multiplied by the Fourier transform of the 
time-reversed series of the Gold code for a satellite 12 
(sampled at the same rate as the pseudorange samples and 

30 stored in memory 216 or 220) . The product is then 
inverse -Fourier transformed to yield the pseudorange time 
sequence . 
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At decisional step 714, a determination is made 
regarding whether a correlation peak exists in each of 
the pseudorange time sequences that is sufficiently above 
a noise-background threshold. According to one 

5 embodiment, a correlation peak is sufficiently above the 
noise -background threshold when the correlation peak is 
at least 6dB above the noise -background threshold. 
However, it will be understood that a correlation peak 
may be sufficiently above the noise -background threshold 

10 when the correlation peak is any suitable level above the 
noise -background threshold. If such a peak does not 
exist for any of the pseudorange time sequences, the 
method follows the No branch from decisional step 714 to 
decisional step 716. 

15 At decisional step 716, a determination is made 

regarding whether a maximum number of combined 
pseudorange samples has been exceeded for the pseudorange 
time sequences without a correlation peak. According to 
various embodiments, the maximum number of combined 

2 0 samples is based on a corresponding amount of sampling 
time for a location estimate that is approximately 1 to 
30 seconds. However, it will be understood that the 
maximum number of combined samples may be any suitable 
number . 

25 If the maximum number of combined pseudorange 

samples has been exceeded, this indicates that the 
corresponding pseudorange signals being received at the 
ranging receiver 14 are too weak to be used. In this 
situation, the method follows the Yes branch from 

30 decisional step 716 to step 717 where a flag is set to 
indicate that the signal is too weak, at which point the 
method comes to an end. However, if the maximum number 
of combined pseudorange samples has not been exceeded, 
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the method follows the No branch from decisional step 716 
and returns to step 70 0 to process more samples to be 
combined with previously compensated and stacked samples 
for pseudorange estimation. 
5 Returning to decisional step 714, if a correlation 

peak that is sufficiently above the noise -background 
threshold does exist for each of the pseudorange time 
sequences, the method follows the Yes branch from 
decisional step 714 to step 718. At step 718, a 

10 pseudorange is determined based on the correlation peak, 
after which the method comes to an end. According to one 
embodiment, a mult i -path correcting centroid calculation 
is applied to determine a substantially exact location of 
the center of each of the peaks based on techniques that 

15 compensate for the distortion of each peak's shape by the 
presence of multi-path signals occurring near the 
correlation peak. In this way, a useable pseudorange may 
be determined based on approximately 0 . 1 seconds worth of 
received signal samples comprising a minimum carrier-to- 

2 0 noise density ratio of approximately 20dB. 

FIGURES 8A and 8B are flow diagrams illustrating a 
method for aligning message data to determine a 
pseudorange in the ranging receiver 14 in a stand-alone 
mode in accordance with one embodiment of the present 

25 invention. The method of FIGURE 8A begins with the 
availability of a specified number of identified 
significant carriers, each having a CFO. According to 
one embodiment, the carriers are identified through the 
direct CFO extraction process described above in 

30 connection with FIGURES 6A and 6B. 

For each of the CFOs, a pointer-step increment, 
g(m), is generated at step 800. This pointer-step 
increment is used to increment a pointer that indexes a 
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table of sine and cosine trigonometric values to be used 
for frequency compensation by complex- vector rotation. 
The magnitude of the increment is proportional to the 
frequency to be compensated. In one embodiment, by using 
5 the moduli and signs of the pointer value, the table can 
be restricted to one quadrant of values. 

At step 8 02, a number of epochs per sub- stack, P, is 
initialized and its counter, p, is set to a value of one. 
According to one embodiment, P is initialized to a value 

10 of four or five. At step 804, a number of sub-stacks per 
analysis, Q, is initialized and its counter, q, is set to 
a value of one . According to one embodiment Q is 
initialized to a value of 60/P, or twelve for the 
embodiment in which P is initialized to a value of five. 

15 At step 806, a number of samples per epoch, K, is 

initialized, and a samples per epoch counter, k, is set 
to a value of one. According to one embodiment, K is 
initialized to a value of 1,023 x 2", with n = 1, 2 or 3 . 
At step 808, a number of carriers, M, is initialized and 

20 a signal sample counter, i, is set to a value of one. 
According to one embodiment, M is initialized to a value 
of four, five, or six, based on the specified number of 
significant carriers available, as described above. At 
step 810, a carrier counter, m, is set to a value of one. 

25 At step 812, the ranging receiver 14 gets a signal 

sample. At step 814, the ranging receiver 14 gets the m^^ 
CFO compensation term for the i*^*^ sample. In one 

embodiment, the i*^^ pointer value for the m"^^ CFO is n = 
modulo [g(m) i] , where the modulus is equal to the number 

3 0 of entries in the trigonometric table that represent the 
interval 0 to 90 degrees of rotation. The sine and 
cosine values at this pointer value are fetched from the 
table to form the complex rotational compensation vector 
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w = cosine [n] + j sine [n] , where j represents the square 
root of -1. At step 816, the ranging receiver 14 complex 
multiplies the sample by the CFG compensation term. At 
step 818, the ranging receiver 14 accumulates the result 
5 of the complex multiplication at a location identified by 
k,q (i.e., position k in sub-stack q) in an m'^*^ array. At 
step 820, the ranging receiver 14 increments m. 

At decisional step 822, a determination is made 
regarding whether or not m is greater than M. If m is 

10 not greater than M, the method follows the No branch from 
decisional step 822 and returns to step 814. However, if 
m is greater than M, the method follows the Yes branch 
from decisional step 822 to step 824. At step 824, the 
ranging receiver 14 increments both k and i. 

15 At decisional step 826, a determination is made 

regarding whether or not k is greater than K. If k is 
not greater than K, the method follows the No branch from 
decisional step 826 and returns to step 810 where m is 
set to a value of one before another signal sample is 

20 collected. However, if k is greater than K, the method 
follows the Yes branch from decisional step 82 6 to step 
828. At step 828, k is set to a value of one. At step 
83 0, the ranging receiver 14 increments p. 

At decisional step 832, a determination is made 

25 regarding whether or not p is greater than P. If p is 
not greater than P, the method follows the No branch from 
decisional step 832 and returns to step 810 where m is 
set to a value of one before another signal sample is 
collected. However, if p is greater than P, the method 

3 0 follows the Yes branch from decisional step 832 to step 
834. At step 834, p is set to a value of one. At step 
836, the ranging receiver 14 increments q. 
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At decisional step 838, a determination is made 
regarding whether or not q is greater than Q. If q is 
not greater than Q, the method follows the No branch from 
decisional step 838 and returns to step 810 where ra is 
5 set to a value of one before another signal sample is 
collected. However, if q is greater than Q, the method 
follows the Yes branch from decisional step 83 8 and comes 
to an end. At this point, sub-stacks have been generated 
for use in determining message data alignment, as 

10 described in more detail below and in connection with 
FIGURES 9, 10 and 11. 

The method of FIGURE 8B begins at step 844 where the 
ranging receiver 14 sets k, p, q and a sub-stacks per bit 
counter, r, all to a value of one and initializes a 

15 number of sub- stacks per data bit, R. According to one 
embodiment, R is initialized to a value of 20/P. At step 
846, the ranging receiver 14 sets m to a value of one. 
At step 848, the ranging receiver 14 adds the value at 
the location identified by k,q of the m'^*' array to the 

20 value at the location identified by k,r of the m^*" array. 
At step 850, the ranging receiver 14 adds the value at 
the location identified by k,q+R of the m'^'^ array to the 
value at the location identified by k,r+R of the m'^*' 
array. At step 852, the ranging receiver 14 increments 

25 m. 

At decisional step 854, a determination is made 
regarding whether or not m is greater than M. If m is 
not greater than M, the method follows the No branch from 
decisional step 854 and returns to step 848 to add the 
3 0 values at the incremented m^^ array. However, if m is 
greater than M, the method follows the Yes branch from 
decisional step 854 to step 856. At step 856, the 
ranging receiver 14 increments k. 
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At decisional step 858, a determination is made 
regarding whether or not k is greater than K. If k is 
not greater than K, the method follows the No branch from 
decisional step 858 and returns to step 846 where m is 
5 set to a value of one. However, if k is greater than K, 
the method follows the Yes branch from decisional step 
858 to step 860. At step 860, the ranging receiver 14 
sets k to a value of one. At step 862, the ranging 
receiver 14 increments p. 

10 At decisional step 864, a determination is made 

regarding whether or not p is greater than P. If p is 
not greater than P, the method follows the No branch from 
decisional step 864 and returns to step 846 where m is 
set to a value of one. However, if p is greater than P, 

15 the method follows the Yes branch from decisional step 
864 to step 866. At step 866, the ranging receiver 14 
sets p to a value of one. At step 868, the ranging 
receiver 14 increments r. 

At decisional step 870, a determination is made 

20 regarding whether or not r is greater than R. If r is 
not greater than R, the method follows the No branch from 
decisional step 870 and returns to step 846 where m is 
set to a value of one. However, if r is greater than R, 
the method follows the Yes branch from decisional step 

25 870 to step 872. At step 872, the ranging receiver 14 
sets both r and k to a value of one. At step 874, the 
ranging receiver 14 sets m to a value of one. 

At step 876, the ranging receiver 14 adds the value 
at the location identified by k,r of the m^^ array to the 

30 value at the location identified by k,r of an m^*^ sum 
stack and to the value at the location identified by 
k,r+R of an m*^^ difference stack. 
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At step 878, the ranging receiver 14 adds the value 
at the location identified by k,r+R of the m'^'^ array to 
the value at the location identified by k,r of the m'^'^ sum 
stack. At step 880, the ranging receiver 14 subtracts 
5 the value at the location identified by k,r+R of the m''*' 
array from the value at the location identified by k,r+R 
of the m''^ difference stack. At step 882, the ranging 
receiver 14 increments m. 

At decisional step 884, a determination is made 

10 regarding whether or not m is greater than M. If m is 
not greater than M, the method follows the No branch from 
decisional step 884 and returns to step 874 where m is 
set to a value of one. However, if m is greater than M, 
the method follows the Yes branch from decisional step 

15 884 to step 886. At step 886, the ranging receiver 14 
increments k. 

At decisional step 888, a determination is made 
regarding whether or not k is greater than K. If k is 
not greater than K, the method follows the No branch from 
20 decisional step 888 and returns to step 874 where m is 
set to a value 1. However, if k is greater than K, the 
method follows the Yes branch from decisional step 888 to 
step 890. At step 890, the ranging receiver 14 sets k to 
a value of one . 

25 At step 892, the ranging receiver 14 correlates each 

sum stack (evaluating the same-signed data bits) and 
difference stack (evaluating the • different-signed data 
bits) with the corresponding Gold codes to obtain 
pseudorange peaks and stores the peak values and their 

30 polarities in an output array at locations (m,r). At 
step 894, the ranging receiver 14 increments r. 

At decisional step 896, a determination is made 
regarding whether or not r is greater than R. If r is 
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not greater than R, the method follows the No branch from 
decisional step 896 and returns to step 874 where m is 
set to a value of one. However, if r is greater than R, 
the method follows the Yes branch from decisional step 
5 896 to step 898. At step 898, the ranging receiver 14 
evaluates the results of the pseudorange correlations 
stored in step 892 in order to determine the quality of 
the message data alignment obtained during the current 
iteration. Based on the maxima of the results, the 

10 ranging receiver 14 determines the alignment, and based 
on the sign of the maxima, the ranging receiver 14 
determines the polarity of the corresponding message data 
bit . At this point , stacks have been generated and used 
to determine coarse message data alignment, as described 

15 in more detail below in connection with FIGURES 9A and 
. 9B. 

FIGURES 9A and 9B are block diagrams illustrating a 
system 900 for generating accumulated sub- stacks and 
stacks of CFO- compensated samples and for finding which 

2 0 of the alignments of those stacks with the message data 
bits produces the largest pseudorange estimates in 
accordance with one embodiment of the present invention. 
By determining alignment for the satellite message data, 
the ranging receiver 14 may remove the satellite 

25 navigation message data modulation from the received 
signal by inverting the signal (multiplying the received 
signal by negative one (-1) ) where the bits of the 
satellite navigation message data modulation result in 
the Gold code modulation being inverted. Multiple, 

30 synchronous repetitions of the signals (i.e., code-epoch 
duration sets of signal samples) may then be added 
together to improve the signal-to-noise ratio. In this 
way, processing gain may be increased as compared to a 
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signal that is squared in order to remove the satellite 
message data. 

The system 900 comprises a data stream 902 that 
represents the signal being received at the ranging 
receiver 14 from the satellites 12. In the illustrated 
embodiment, each message data bit comprises 20 code 
epochs. For the illustrated example, a portion of a 
first message data bit 902a is followed by a second 
message data bit 902b, a third message data bit 902c and 
a fourth message data bit 902d. The data stream 902 
comprises a plurality of epochs 912 that are received at 
the ranging receiver 14 from the satellites 12. In the 
illustrated embodiment, for GPS signals, there are 1,023 
Gold-code chips in an epoch 912, which may be sampled at 
2,048 or 4,096 samples per code epoch or other suitable 
sampling rate. Thus, each epoch 912 may include 2,048 
signal samples for the first sampling rate. 

TO determine the satellite message data alignment in 
accordance with one embodiment of the present invention, 
the ranging receiver 14 combines a plurality of code- 
epoch duration sets of samples 912 into a plurality of 
sub-stacks 914. According to the illustrated embodiment, 
each sub-stack 914 comprises five code-epoch duration 
sample sets 912. The ranging receiver 14 then further 
, combines the sub- stacks 914 into a plurality of stacks 
916. According to the illustrated embodiment, each stack 
916 comprises four sub-stacks 914. In an alternative 
embodiment, each sub- stack 914 could comprise four code 
epoch duration sample sets 912 and each stack 916 could 
0 comprise five sub-stacks 914. It will be understood that 
other suitable numbers of code-epochs duration sample 
sets 912 may be included in each sub-stack 914 and other 
suitable numbers of sub-stacks 914 may be included in 



DALOl: 756059.1 



ATTORNEY'S DOCKET 
068796 - 0103 



56 



PATENT APPLICATION 



each stack 916 without departing from the scope of the 
present invention. 

The ranging receiver 14 then identifies pairs of 
stacks 916 such that each pair comprises adjacent code- 
5 epoch duration sample sets 912 of the data stream 902. 
For example, the pair 916a comprises stack SI, which 
includes the first 20 code-epoch duration sample sets 912 
of the data stream 902, and stack S7, which includes the 
second 20 code-epoch duration sample sets 912 of the data 

10 stream 902. For each such pair of stacks 916, the 
ranging receiver 14 forms both the sum and difference 
values in each of the stacks 916. The message data 
alignment may then be determined based on the largest 
pseudorange peak values for each sum and difference pair. 

15 In addition, the polarities of the largest pseudorange 
peak values indicate the signs for the corresponding 
message data bit pairs. Knowledge of whether the sum or 
difference provided the larger pseudorange peak value may 
be used to determine the order of the data bit-pair 

20 {e.g., "01" or "10" for difference and "00" or "11" for 
sum) . 

FIGURE 9B is a block diagram illustrating a portion 
of the system 900 in greater detail in accordance with 
one embodiment of the present invention. The illustrated 

25 system 900 comprises a plurality of code-epoch duration 
sample sets 912a to 912e that comprise signal samples 
that are representative of the Gold code chips for the 
signal received from a satellite 12 at the ranging 
receiver 14 . To determine coarse satellite message data 

3 0 alignment, the ranging receiver 14 begins at an arbitrary 
time to accumulate samples of the signal for a satellite 
12 in accumulators 918. 
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As used herein, "coarse" alignment means alignment 
to within 2.5 epochs 912, which arises from the choice of 
5 epochs 912 per sub-stack 914. According to one 
embodiment, the number of accumulators 918 corresponds to 
5 the number of samples in an epoch 912. Thus, for the 
embodiment in which each epoch 912 comprises 1,023 chips 
and the signal -sampling rate is 2,048 samples per code 
epoch, the system 900 comprises 2,048 accumulators 918. 

As each sample of the signal is received at the 

10 ranging receiver 14, the signal is accumulated in an 
accumulator 918. For example, if the sampling rate was a 
sample per chip and synchronous with its position in the 
Gold code epoch and if the ranging receiver 14 begins to 
accumulate samples of the Gold code with the 733^"^ chip of 

15 the Gold code, the sample representing the 73 3^*^ chip of 
the Gold code would be accumulated in the Ai accumulator 
918, the sample representative of the 734*^'* chip of the 
Gold code would be accumulated in the A2 accumulator 918, 
and so on, with the sample representing the 732"** chip of 

20 the Gold code accumulated in the An accumulator 918. 

After the samples of the first code-epoch duration 
sample set 912, which comprises an entire cycle of the 
Gold code, are accumulated in the accumulators 918, the 
samples of the next code-epoch duration sample set 912, 

25 or cycle of samples of the Gold code signal, is similarly 
accumulated such that each accumulator 918 accumulates 
the sample representing the same chip of the Gold code. 
Thus, while within the same navigation message data bit, 
the samples representing the first chip of each epoch 912 

3 0 are accumulated with the same algebraic sign, and so on 
for all the samples in the epochs 912 representing the 
current message data bit . When the sign of the next 
message bit changes, the sign of all the signal samples 
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representing the new message data bit are inverted and, 
when accumulated, diminish the previously accumulated 
sums . 

According to the illustrated embodiment, each 
accumulator 918 accumulates samples for five code-epoch 
duration sample sets 912. However, it will be understood 
that the accumulators 918 may accumulate samples for any 
suitable number of epochs 912 without departing from the 
scope of the present invention. 

After five code-epoch duration sample sets 912 are 
accumulated in the accumulators 918, the resulting values 
in the accumulators 918 are saved into a sub-stack 914, 
such as the sub-stacks 914 labeled A through M in the 
illustrated embodiment. The accumulators 918 may then be 
cleared and begin accumulating samples for the duration 
of the next five epochs 912 to be combined into the next 
sub- stack 914. 

For the embodiment in which the accumulators 918 
accumulate chips for the duration of five epochs 912, 
each sub- stack 914 comprises a combination of five code- 
epoch duration sample sets 912. However, it will be 
understood that the sub- stacks 914 may comprise any 
suitable number of code-epoch duration sample sets 912 
based on the number of code-epoch durations accumulated 
in the accumulators 918 without departing from the scope 
of the present invention. 

The sub- stacks 914 are further combined into stacks 
916. According to one embodiment, each stack 916 
comprises four sub-stacks 914. Thus, for the embodiment 
in which each sub-stack 914 comprises five code-epoch 
duration sample sets 912, each stack 916 comprises 20 
code-epoch duration sample sets 912. This corresponds to 
the length of a bit of satellite navigation message data. 
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which is 20 epochs 912. However, it will be understood 
that the stacks 916 may comprise any suitable number of 
sub- stacks 914 without departing from the scope of the 
present invention. 
5 The stacks 916 are grouped into pairs of stacks 916 

representing adjacent data bits. Thus, for example, the 
pair 916a comprises stack SI (comprising sub-stacks A, B, 
C and D) and stack S7 (comprising sub-stacks E, F, G and 
H) . The pair 916b comprises stack S2 (comprising sub- 

10 stacks B, C, D and E) and stack S8 (comprising sub-stacks 
F, G, H and I) . For the illustrated embodiment based on 
sub-stacks 914 of A through M, the sub-stacks 914 may be 
combined into twelve stacks 916, or six pairs of stacks 
916a, 916b, 916c, 916d, 916e and 916f. 

15 To determine a coarse alignment for the satellite 

message data, the stacks 916 in each pair are summed and 
differenced and the sum and difference arrays are 
evaluated for pseudorange peaks. Thus, for example, the 
values of Stack SI and Stack S7 are added together and 

2 0 subtracted from each other, while the values of Stack S2 

and Stack SB are added together and subtracted from each 
other, and so on for each of the six pairs of stacks 
916a, 916b, 916c, 916d, 916e and 916f. 

FIGURE 10 is a block diagram illustrating a system 
25 920 for refining satellite message data alignment in the 
system 10 in accordance with one embodiment of the 
present invention. The refining process begins with the 
two stacks 916 derived in the coarse alignment method 
described above in connection with FIGURE 9. 

3 0 For coarse alignment to have occurred, the two 

message data bits represented by the stacks producing the 
match, such as SI and S7, are either of opposite logical 
value (e.g., ..x,0,l,x.. or ..x,l,0,x-.) or are bounded 
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by bits of opposite logical value (e.g., ..0,1,1,0.. or 
..1,0,0,1..). These two stacks 916 are sub-divided into 
partial stacks 922a and 922b and refining stacks 924 V, 
W, X and Y. 

5 Two refining stacks 924 U and Z may be included to 

provide an enlarged refining shift range with higher 
resolution than was used in the coarse alignment method 
described above in connection with FIGURE 9. This 
refining shift range may allow the stacks 916 to be 

10 shifted earlier by an interval represented by two epochs 
912 and later by the same amount of time for a more 
refined alignment. According to one embodiment, the 
refining stacks 924 are formed from the signal samples 
912, which have been saved previously during the coarse 

15 alignment stacking process, as described further below. 

Given that the choice of a sub- stack duration of 
five epochs 912 in the coarse alignment process described 
above provided an alignment of within 5/2, or 2,5, epochs 
912, the choice of two refining stacks 924 in this 

20 illustration tests the alignment further by two refining 
stacks 924 and moves the alignment within 1 epoch 912, 
which is referred to herein as a "refined" alignment. As 
compared to the system 900, this system 920 reduces 
computational load for reaching a closer alignment by 

25 adding or subtracting only small ref ining-stacks 924 to 
the ends of the coarse alignment stacks 916 derived in 
the system 900. 

According to one embodiment, each partial stack 922 
comprises sixteen code epochs' worth of CFO-compensated 

30 signal samples and each refining stack 924 comprises two 
code epochs' worth. Thus, according to the embodiment in 
which each stack 916 comprises 20 code epochs' worth of 
samples, a stack 916 comprises a partial stack 922 and 
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two refining stacks 924. For the illustrated embodiment, 
a refined alignment stack 916 may comprise the partial 
stack 922a with the refining stacks 924 of U and V, V and 
W, or W and X, i.e., allowing a total of three possible 
5 shifts spanning a total of four epochs 912. The 
corresponding stack 916 for a pair of stacks 916 may 
comprise the partial stack 922b with the refining stacks 
924 of W and X, X and Y, or Y and Z. 

According to one embodiment, the ranging receiver 14 

10 synchronously accumulates samples of the signal from 
satellites 12 in the accumulators 918 as described above 
in connection with FIGURE 9B. However, in the event that 
refined alignment will be performed, pairs of code -epoch 
duration sample sets formed by the system 900 during the 

15 coarse alignment process are temporarily stored for later 
use in the refined alignment process. In this case, when 
the coarse alignment has been determined, additional 
refining stacks 924 U, V, W, X, Y and Z are formed out of 
the corresponding sample-set pairs saved. 

20 For this embodiment, the number of epoch summations 

to be performed to reach a particular alignment quality 
is reduced over the number required if the system 900 is 
used. First, summing each of the six refining stacks 924 
results in six sums. Second, the aligning stacks 916 may 

25 be shifted by adding or subtracting the appropriate 
refining stacks 924 to shift the alignment in the desired 
direction . 

For example, if the desired shift is to the left 
(earlier) , then the shifted SI stack 916 would include SI 
30 + U - W, and the shifted 57 stack 916 would include 87 + 
W - Y. Once these shifted stacks 916 have been formed, 
the shifted stacks 916 are evaluated as described in 
relation to the system 900 above. If the desired shift 
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is to the right (later) , then the shifted SI stack 916 
would include 81 + X - V, and the shifted 37 stack 916 
would include 37 + Z - X, and the outcome would again be 
evaluated as described in relation to the system 900 
5 above. Thus, only two additional single-epoch duration 
adds and subtracts per stack pair, over the computations 
carried out in the coarse alignment process using the 
system 900, allows the alignment to be refined to within 
about one epoch 912. 

10 This process of alignment refinement by smaller 

shifts may be extended to achieve any suitable degree of 
alignment refinement, but alignment to within less than 
an epoch's duration yields very little additional 
processing gain. 

15 Once a coarse or, if desired, a refined alignment 

has been achieved, the alignment may be used to decode 
the continuing satellite message data, two message data 
bits at a time, by carrying out a simpler, 20 plus 20 
code -epoch stacking routine for each of the subsequent 

2 0 bit pairs. The magnitude and polarity of the sum or 

difference determines the values of the respective 
message data bits. 

FIGURE 11 is a flow diagram illustrating a method 
for aligning message data in the ranging receiver 14 in 
25 accordance with one embodiment of the present invention. 
The method begins at step 1100 where the ranging receiver 
14 generates sub-stacks 914 based on epochs 912 
accumulated in the accumulators 918. At step 1102, the 
ranging receiver 14 combines the sub- stacks 914 into 

3 0 stacks 916 which are grouped together in pairs. At step 

1104, the ranging receiver 14 calculates a sum and 
difference for each pair of stacks 916 . 
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At Step 1106, the ranging receiver 14 determines the 
maxima of the pseudorange peaks, the peaks being based on 
correlations of the stacks 916 with the satellite Gold 
codes. According to one embodiment, this is accomplished 
5 by searching for the six largest values in the correlates 
of the sum and difference stacks, along with the offsets 
at which the maxima occur. After the maxima are found, a 
satellite ID is determined from the Gold code yielding 
the peak and the pseudorange is determined by the 
10 centroid of the correlation peak. In this way, the 
satellite message data may be aligned to within 2.5 
epochs 912 . 

At decisional step 1108, a determination is made 
regarding whether or not the alignment is to be refined. 
15 If the alignment is not to be refined, the method follows 
the No branch from decisional step 1108 and comes to an 
end. However, if the alignment is to be refined, the 
method follows the Yes branch from decisional step 1108 
to step 1109. 

20 At step 1109, the ranging receiver 14 identifies 

partial stacks 922 and refining stacks 924 which may be 
combined into stacks 916. At step 1110, the ranging 
receiver 14 determines sums for the refining stacks 924. 
At step 1112, the ranging receiver 14 shifts the stacks 

25 916 based on the refining stacks 924, as described in 
more detail above in connection with FIGURE 10. At step 
1114, the ranging receiver 14 determines a refined 
alignment based on the shifted stacks 916, at which point 
the method comes to an end. In this way, the satellite 

30 message data may be aligned to within 0.5 epochs 912. 

FIGURE 12 is a flow diagram illustrating a method 
for predicting satellite message data received in 
positioning signals at the ranging receiver 14 in 
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accordance with one embodiment of the present invention. 
The method begins at decisional step 1200 where a 
determination is made regarding the status of a 
superf rame in the ranging receiver 14 . 
5 If no superframe is present in the ranging receiver 

14 or if the superframe is expired, the method follows 
the Expired branch from decisional step 1200 to step 
1202. The superframe may be considered expired if the 
superframe is greater than 96 hours old. However, it 

10 will be understood that the superframe may be considered 
expired based on another suitable amount of time or other 
suitable criteria without departing from the scope of the 
present invention. At step 1202, the ranging receiver 14 
downloads the superframe from the satellites 12 or from 

15 the geolocation processor 24 . 

If the superframe in the ranging receiver 14 is not 
current, the method follows the Not Current branch from 
decisional step 1200 to step 1204. The superframe may be 
considered not current if the superframe is between 4 and 

20 96 hours old. However, it will be understood that the 
superframe may be considered not current based on another 
suitable amount of time or other suitable criteria 
without departing from the scope of the present 
invention. At step 1204, the ranging receiver 14 

25 downloads an update to the superframe from the satellites 
12 or from the geolocation processor 24. 

Any data obtained from the geolocation processor 24 
in steps 1202 or 1204 may be obtained over the voice 
channel using the modems 40 or via an alternate pathway 

3 0 using, for example, a control channel for the wireless 
carrier network 2 0 or other suitable data pathway. 

Returning to decisional step 1200, if the superframe 
in the ranging receiver 14 is current, the method follows 
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the current branch from decisional step 1200 to 
decisional step 1206. Also, after a superframe has been 
downloaded in step 1202 or an update has been downloaded 
in step 1204, the method continues to decisional step 
1206. The superframe may be considered current if the 
superframe is less than 4 hours old. However, it will be 
understood that the superframe may be considered current 
based on another suitable amount of time or other 
suitable criteria without departing from the scope of the 
present invention. 

At decisional step 1206, a determination is made 
regarding whether or not the real-time clock 232 has the 
current time. If the real-time clock 232 does not have 
the current time, the method follows the No branch from 
decisional step 1206 to step 1208. At step 1208, the 
real-time clock 232 for the ranging receiver 14 is 
calibrated, or synchronized, to GPS time based on a 
remote clock or TOD information in the positioning 
signals. Once calibrated, the real-time clock 232 keeps 
track of GPS time. 

At step 1210, the ranging receiver 14 sets up a 
common table that is operable to store data that is 
common to the superframe for each satellite in the 
constellation of satellites. According to one 

embodiment, approximately 60% of the superframe is common 
to each satellite in the constellation. 

At step 1212, for each of the satellites, the 
ranging receiver 14 sets up a unique table that is 
operable to store data that is unique to the superframe 
, of the corresponding satellite. For example, this unique 
data may comprise ephemeris data, telemetry data, 
tropospheric data, ionospheric propagation models, and 
other suitable data relating to the satellite. At this 



DALOl: 756059 .1 



ATTORNEY'S DOCKET 
068796.0103 



66 



PATENT APPLICATION 



point, along with time and date information from the 
real-time clock 232 at the time it is required, the 
message predictor has been set up for subsequent use by 
the ranging receiver 14 in predicting satellite message 
5 data for a specified satellite 12 from which positioning 
signals are received. From step 1212, the method 
continues to step 1214. 

Returning to decisional step 1206, if the real-time 
clock 23 2 has the current time, the method follows the 

10 Yes branch from decisional step 12 06 to decisional step 
1215. At decisional step 1215, a determination is made 
regarding whether the prediction tables, which include 
the common table and the unique tables, are current. 
Each of the prediction tables may comprise any suitable 

15 data store. If the prediction tables are not current, 
the method follows the No branch from decisional step 
1215 and proceeds to step 1210 where the ranging receiver 
14 sets up a common table. 

However, if the prediction tables are current, the 

20 method follows the Yes branch from decisional step 1215 
to step 1214. At step 1214, the ranging receiver 14 
begins to predict a message fragment. At step 1216, the 
ranging receiver 14 identifies the satellite 12, and the 
time interval, corresponding to the message fragment to 

25 be predicted. At step 1218, the ranging receiver 14 
generates a message fragment table for predicting the 
message fragment . 

At step 1220, the ranging receiver 14 adds data from 
the common table to the message fragment table. At step 

30 1222, the ranging receiver 14 adds data from the unique 
table for the satellite 12 identified in step 1216 to the 
message fragment table. At step 1224, the ranging 
receiver 14 adds time -dependent data, such as bit 
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transition timing information, to the message fragment 
table based on the current GPS time obtained from the 
real-time clock 232. 

At step 1226, the ranging receiver 14 provides the 
5 predicted message fragment, which comprises the data for 
predicting the satellite message data, to the requester. 
Thus, in this way, the ranging receiver 14 may predict 
the satellite message data in the positioning signals so 
that the message data may be removed after alignment, 

10 which may be performed as described above in connection 
with FIGURES 9, 10 and 11. 

FIGURES 13A, 13B and 13C are flow diagrams 
illustrating a method for determining a pseudorange 
estimate in the ranging receiver 14 in a geometric mode 

15 in accordance with one embodiment of the present 
invention. The method provides for determining a 

pseudorange based on the availability at the ranging 
receiver 14 of: (i) accurately known real time (with 
reference to the satellite constellation and within about 

20 a millisecond) ; (ii) an estimate of the current location 
of the ranging receiver 14 (within about a mile) ; 
(iii) an estimate of the frequency bias of the ranging 
receiver 14; (iv) recent (within about 24 hours) 
ephemeris data; and (v) recent (within about 96 hours) 

25 constellation superframe data. 

For example, real time data may be received from a 
first set of satellites, ephemeris data may be received 
from a second set of satellites, and superframe data may 
be received from a third set of satellites. The sets of 

30 satellites, which may each comprise one or more 
satellites, may comprise none, some or all of the same 
satellites in the satellite constellation. According to 
one embodiment, the first set of satellites comprises one 

DALOl: 756059.1 



ATTORNEY'S DOCKET 
068796 . 0103 



68 



PATENT APPLICATION 



of the in-view satellites 12, the second set of 
satellites comprises at least four of the in-view 
satellites 12, and the third set of satellites comprises 
one of the satellites in the constellation. In an 
5 alternative embodiment, the application server 64 and/or 
the workstation 46 may serve as a satellite information 
server from which ephemeris data and/or superframe data 
may be received. 

In FIGURE 13A, the method begins at step 130 0 where 

10 the ranging receiver 14 determines a set of Doppler 
frequency shift, Doppler-rate and satellite ID data for a 
selected set of satellites from calculations based on the 
current geometric data referred to above. This set of 
data may be selected from among the satellites 12 

15 potentially viewable from the location of the ranging 
receiver 14 at that time. This selection may be based on 
suitable criteria, such as low geometric dilution of 
precision, low elevation angle, or any suitable criteria. 
At step 1302, for each satellite 12 in the set, the 

20 ranging receiver 14 sets up a table of message data bits 
that will be transmitted during the time that pseudorange 
samples will be processed. 

At step 13 04, the ranging receiver 14 computes the 
CFG and CFO-rate for each of the satellites 12. The CFG 

25 is determined by adding the frequency bias of the ranging 
receiver 14 to the computed Doppler frequency shift. The 
ranging receiver 14 then assigns the Doppler rate to the 
CFO-rate for each of the satellites 12. At step 1306, 
the ranging receiver 14 stores the CFGs and the 

30 corresponding CFO-rates in a table. 

At step 1308, the ranging receiver 14 initializes a 
maximum number of samples, I, and sets a sample counter, 
i, to one. The maximum number of samples may be chosen 
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to limit the maximum integration time to suit the 
application, such as 10 to 30 seconds, but other values 
outside this range are also possible. At step 1310, a 
number of epochs, P, to be averaged is initialized. The 
5 number of epochs (which corresponds to the number of 
samples) to be averaged during one loop-through may be 
selected according to various criteria, such as the total 
processing gain to be provided per loop- through, the 
amount of memory available for table or intermediate 

10 result storage and any other suitable criteria. At step 
1312, a number of carriers, M, is initialized. 

At step 1314, a number of samples per epoch, K, is 
initialized. At step 1316, the ranging receiver 14 sets 
the epoch sample number, k, to a value of one. At step 

15 1318, the ranging receiver 14 sets the carrier counter, 
m, to a value of one. 

At step 1320, the ranging receiver 14 gets a signal 
sample for processing. At step 1322, the ranging 
receiver 14 gets the m*^*' CFO compensation term for the i*^*' 

20 sample. At step 1324, the ranging receiver 14 complex 
multiplies the sample by the CFO compensation term. At 
step 1326, the ranging receiver 14 scalar multiplies the 
CFO-compensated sample by its corresponding data bit 
polarity, which may be determined in accordance with the 

25 method of FIGURE 13B. 

At step 1328, the ranging receiver 14 adds the 
result of the multiplication in step 1326 to a location 
identified by k in an m*^^ array. At step 1330, the 
ranging receiver 14 increments m. At decisional step 

30 1332, a determination is made regarding whether or not m 
is greater than M. If m is not greater than M, the 
method follows the No branch from decisional step 1332 
and returns to step 1322. However, if m is greater than 
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M, the method follows the Yes branch from decisional step 
1332 to step 1334. At step 1334, the ranging receiver 14 
increments both k and i. 

At decisional step 1336, a determination is made 
5 regarding whether or not k is greater than K. If k is 
not greater than K, the method follows the No branch from 
decisional step 1336 and returns to step 1318 where m is 
set to a value of one. However, if k is greater than K, 
the method follows the Yes branch from decisional step 

10 1336 to step 1338. At step 1338, the ranging receiver 14 
sets k to a value of one. At step 134 0, the ranging 
receiver 14 increments p. 

At decisional step 1342, a determination is made 
regarding whether or not p is greater than P. If p is 

15 not greater than P, the method follows the No branch from 
decisional step 1342 and returns to step 1318 where m is 
set to a value of one. However, if p is greater than P, 
the method follows the Yes branch from decisional step 
1342 to step 1344. At step 1344, the ranging receiver 14 

20 sets m to a value of one. 

At step 1346, the ranging receiver 14 correlates an 
m*^*^ pseudorange stack with a Gold code for an m^^ satellite 
ID. At step 1348, the ranging receiver 14 compensates 
for time dilation by correlation peak shifting. According 

25 to one embodiment, time-dilation compensation may be 
accomplished by vector rotation in the frequency domain 
during the correlation by FFT means during step 1346. 
According to another embodiment, time-dilation 
compensation may be accomplished by suitable sub-stack 

30 sample shifting prior to final accumulation and 
conversion to pseudorange by correlation. In yet another 
embodiment, the time-dilation compensation may be 
performed in accordance with the time-shifting method of 
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FIGURE 13C. At step 1350, the pseudorange stack is added 
to an m''*' extended output stack and the m"'^ pseudorange is 
evaluated for an adequate peak. At step 1352, the 
ranging receiver 14 increments m. 
5 At decisional step 1354, a determination is made 

regarding whether or not m is greater than M. If m is 
not greater than M, the method follows the No branch from 
decisional step 1354 and returns to step 1346 where the 
m*^^ pseudorange stack is correlated with the Gold code of 

10 the m*^^ satellite ID for the incremented m. However, if m 
is greater than M, the method follows the Yes branch from 
decisional step 1354 to decisional step 1356. 

At decisional step 13 56, a determination is made 
regarding whether any peak is not adequate. If each peak 

15 is adequate, then the minimum number, M, of pseudoranges 
has been detected successfully and the method follows the 
No branch from decisional step 13 56 and comes to an end. 
In this way, a pseudorange may be determined in the range 
of 1 to 30 seconds based on received signals comprising a 

20 minimum carrier-to-noise density ratio in the range of 
5dB to lOdB. However, if any peak is found to be not 
adequate, the method follows the Yes branch from 
decisional step 1356 to decisional step 1358. 

At decisional step 1358, a determination is made 

25 regarding whether or not a maximum number of samples, I, 
has been exceeded. If the maximum number of samples has 
been exceeded, the method follows the Yes branch from 
decisional step 1358 to step 1360. At step 1360, a flag 
is set to indicate that the signal is too weak for the 

3 0 minimum number, M, of pseudoranges to be detected 
successfully, at which point the method comes to an end. 

Returning to decisional step 1358, if the maximum 
number of samples has not been exceeded, the method 
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follows the No branch from decisional step 1358 to step 
1362. At step 1362, the ranging receiver 14 clears the 
input sample stacks. At this point, the method returns 
to step 1316 where k is set to a value of one. 
5 FIGURE 13B is a flow diagram illustrating a method 

for removing message data modulation in the method of 
FIGURE 13A in accordance with one embodiment of the 
present invention. The method begins at step 13 64 where, 
based on geometric prediction, the ranging receiver 14 

10 obtains satellite IDs, Doppler frequency shift, and 
Doppler rates. At step 13 66, the ranging receiver 14 
stores a message bit transition offset, T{m), for each of 
the M satellites 12 based on the results of the method 
for predicting satellite message data described in 

15 FIGURE 12 . These offsets also include the signal transit 
time to the ranging receiver 14 for each satellite 12. 

At step 1368, the ranging receiver 14 stores a 
specified number of message data bits for each satellite 
12 based on the results of the method for predicting 

20 satellite message data described in FIGURE 12. According 
to one embodiment, the specified number of message data 
bits is equal to the result of dividing (a) the maximum 
number of samples by (b) 20 times the product of the 
sampling rate (samples per chip) and the chips per epoch. 

25 However, it will be understood that any suitable number 
of message data bits may be stored without departing from 
the scope of the present invention. Also, according to 
one embodiment, message data bits of logical (1) and (0) 
may be stored as values of +1 and -1, respectively, or as 

30 values of -1 and +1, respectively. However, it will be 
understood that the message data bits may be stored in 
any suitable manner without departing from the scope of 
the present invention. 
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At step 13 70, the ranging receiver 14 multiplies the 
CFO- compensated samples by the corresponding stored 
message data bits. According to one embodiment, for the 
m*^'^ satellite, the value of the pointer to the message 
5 data bit may be determined by the ranging receiver 14 
based on the following formula: 

IntegerPart [{T(m) +i}/K] +1, 
where T (m) is the message bit transition offset for 
the message bits from the m*^*^ satellite 12, i is the 

10 sample number and K is the number of samples per epoch. 

FIGURE 13C is a flow diagram illustrating a method 
for compensating for time dilation in the method of 
FIGURE 13A in accordance with one embodiment of the 
present invention. In this method, a time offset, t, 

15 from the beginning of the analysis to the beginning of 
the current sub- stack has a value of i-PK when the 
routine is entered from step 1348 of FIGURE 13A. The m"'^ 
Doppler at the beginning of this sub- stack, n (m) , is 
equal to d(m)+t r(m}, where d(m) is the m*^^ Doppler at the 

20 beginning of the analysis and r(m) is the m*^^ Doppler rate 
expressed in per-sample terms, rather than per-second 
terms. The average of the m*^^ Doppler, ad (m) , from the 
beginning of the current analysis to the beginning of the 
current sub-stack is equal to d(m)+(t r(m))/2. The 

25 average of the m*^*^ Doppler in the current sub- stack, 
adstack{m), is equal to n(m)+t r(m)/2. The time shift to 
the beginning of this sub-stack, Del (m) , is equal to 
ad(m)Tt/f(c) in second terms and ad(m)t/fc in sample 
terms. The time shift within the current sub-stack, 

30 Dstack(m), is equal to adstack (m) TPK/f (c) in second terms 
and ad(m)PK/fc in sample terms. The total time shift, 
TS(m), is equal to Del (m) +Dstack (m) . 
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In addition, the following variables are defined for 
the purpose of this illustration, without limiting the 
application to only the GPS satellite constellation 
characteristics, as follows: i = sample number; M = 
5 number of carriers; d(m) = value of Doppler frequency 
shift, per instance of m; n(m) = value of Doppler at the 
beginning of the stacking process (i.e., when i=l) + the 
accumulated bias rate up to the beginning of the current 
sub-stack, per instance of m; r(m) = value of Doppler 

10 rate, per instance of m; K = number of samples per epoch; 
P = number of epochs per integration cycle; S = number of 
samples per sub-stack; c = speed of light; f{c) = GPS 
satellite carrier frequency (i.e., 1575.42 MHz); T = 
interval between samples; v = sample index within current 

15 sub- stack; u = sample index within extended output stack; 
TS = total time shift; SS = sub-stack sample; OS = output 
stack sample; FP [] = FractionalPart of; and IP[] 
IntegerPart of . 

The method begins at step 1380 where the ranging 

20 receiver 14 sets the sample index, v, to a value of one. 
At decisional step 13 82, a determination is made 
regarding whether the time shift is positive or negative 
(which depends on, and is opposite to, the sign of the 
Doppler frequency shift) . If the time shift is negative, 

25 the method follows the Negative branch from decisional 
step 1382 to step 1384. At step 1384, the ranging 
receiver 14 adds (FP [TS] ) *SS (v) to OS (m) , (v+IP [TS] -1) . 
At step 1386, the ranging receiver 14 adds (1- 
FP [TS] ) *SS (v) to OS (m) , (v+IP[TS] ) . At step 1388, the 

30 ranging receiver increments v. 

At decisional step 1390, a determination is made 
regarding whether or not v is greater than K. If v is 
not greater than K, the method follows the No branch from 
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decisional step 1390 and returns to step 1384. However, 
if V is greater than K, pseudorange peaks may be 
evaluated and the method follows the Yes branch from 
decisional step 13 90 and returns to step 1350 in 
5 FIGURE 13A. 

Returning to decisional step 1382, if the time shift 
is positive, the method follows the Positive branch from 
decisional step 1382 to step 1392. At step 1392, the 
ranging receiver 14 adds (1-FP [TS] ) *SS (v) to 

10 OS (m) , (v+IP [TS] ) . At step 1394, the ranging receiver 14 
adds FP[TS]*SS(v) to OS (m) , {v+IP [TS] +1) . At step 1396, 
the ranging receiver 14 increments v. 

At decisional step 13 98, a determination is made 
regarding whether or not v is greater than K. If v is 

15 not greater than K, the method follows the No branch from 
decisional step 1398 and returns to step 1392. However, 
if V is greater than K, pseudorange peaks may be 
evaluated and the method follows the Yes branch from 
decisional step 13 98 and returns to step 13 50 in 

2 0 FIGURE 13A. 

FIGURE 14 is a flow diagram illustrating a method 
for determining geolocation in the ranging receiver 14 
with matching assistance in accordance with one 
embodiment of the present invention. The method begins 

25 at step 1400 where the ranging receiver 14 waits in a 
low-power standby mode for a trigger to start the 
pseudorange estimating process. The trigger may comprise 
a request generated by the wireless device 16 
automatically, a request generated by a user of the 

30 wireless device 16 such as by pushing a button 30 on the 
ranging receiver 14 , a request from an external source 
such as the Internet, a geolocation processor 24 or other 
suitable external source, or any other suitable trigger. 
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Once the trigger is received, the location circuitry 
of the ranging receiver 14 is turned on at step 1402. 
Thus, for example, the power management circuitry 23 0 may 
provide full power to each of the components 204, 206, 
5 208, 210, 212, 214, 216, 218, 220 and 222 after receiving 
the trigger. 

At step 1404, the real-time clock 232 for the 
ranging receiver 14 is calibrated, or synchronized, to a 
remote clock across the wireless link or based on TOD 

10 information in the positioning signals. Once calibrated, 
the real-time clock 232 keeps track of real time so that 
the timing of events in the ranging receiver 14 can be 
accurately tracked. 

At step 1406, the ranging receiver 14 clears any 

15 prior signal samples in the signal sample memory 216. At 
step 1408, the ranging receiver 14 begins to collect 
samples for carrier detection by direct CFO extraction 
from the signals received from the satellites 12. At 
step 1410, the CFOs are directly extracted from the 

20 samples, as described in more detail above in connection, 
with FIGURES 6A and 6B. 

At decisional step 1412, a determination is made 
regarding whether at least a specified number of 
significant carriers and their CFOs have been found. 

25 According to one embodiment, the specified number is 
four. However, it will be understood that the specified 
number may comprise any suitable number. If the 

specified number of significant carriers has not been 
found, the method follows the No branch from decisional 

30 step 1412 to decisional step 1414. 

At decisional step 1414, a determination is made 
regarding whether a maximum number of collected CFO 
extraction samples has been exceeded. According to one 
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embodiment, the maximum number of collected samples 
comprises the number of samples that may be collected in 
approximately one second. Thus, the maximum number of 
collected samples may be based on the sampling rate. 
5 However, it will be understood that the maximum number of 
collected samples may be any suitable number. 

If the maximum number of collected CFO extraction 
samples has been exceeded, this indicates that the 
signals being received at the ranging receiver 14 are too 

10 weak to be used. In this situation, the method follows 
the Yes branch from decisional step 1414 and comes to an 
end. However, if the maximum number of collected CFO 
extraction samples has not been exceeded, the method 
follows the No branch from decisional step 1414 and 

15 returns to step 1408 for the collection of more signal 
samples for CFO extraction. 

Returning to decisional step 1412, if the specified 
number of significant carriers has been found, the method 
follows the Yes branch from decisional step 1412 to step 

20 1416. At step 1416, the ranging receiver 14 sends CFO 
results to the geolocation processor 24 through the 
service center 18. According to one embodiment, the CFO 
results comprise the CFO estimates, the real time for the 
CFO estimate measurements and/or other suitable data 

25 related to the CFO estimates. At step 1418, the ranging 
receiver 14 receives satellite information from the 
geolocation processor 24 . The satellite information may 
comprise satellite IDs, the most recent satellite message 
data fragments, and bit transition information for each 

30 of the satellites 12 identified by the geolocation 
processor 24 based on the CFO results and/or any other 
suitable information. 
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At step 1420, the ranging receiver 14 associates the 
CFOs with their respective satellites 12 based on the 
satellite information received from the geolocation 
processor 24. At step 1422, the ranging receiver 14 
5 begins to collect samples for pseudorange estimation from 
the signals received from the satellites 12. At step 
1424, the ranging receiver 14 begins to process the 
pseudorange samples by compensating the samples for the 
CFOs. 

10 At step 1426, the ranging receiver 14 aligns 

satellite message data fragments with the pseudorange 
samples and removes the message data modulation from the 
samples by multiplying out the biphase message data 
modulation, as described in more detail above in 

15 connection with FIGURES 8, 9, 10, 11 and 12. 

At step 1428, the ranging receiver 14 linearly 
combines, or stacks, the samples into a single stack 
result of one-epoch duration for each satellite 12. At 
step 1430, the ranging receiver 14 correlates the single 

20 epoch -duration stack with the Gold codes for the 
satellites 12. 

According to one embodiment, the ranging receiver 14 
correlates the single epoch-duration stack result with 
the Gold codes by Fourier transform. For this 

25 embodiment, the Fourier transform of each CFO- compensated 
stack of samples is multiplied by the Fourier transform 
of the time-reversed series of the Gold code for a 
satellite 12 (sampled at the same rate as the pseudorange 
samples and stored in memory 216 or 220) . The product is 

3 0 then either compensated for time-dilation by complex 
vector multiplication in the frequency domain before it 
is inverse -Fourier transformed to yield the time-dilation 
compensated pseudorange time sequence, or the product is 
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inverse-Fourier transformed to yield the un-compensated 
pseudorange time sequence, which is compensated for time 
dilation by element -shifting, as described in FIGURE 13C, 
at step 1432. 

5 According to an alternative embodiment, the ranging 

receiver 14 correlates the single epoch- duration stack 
result with the respective Gold code in hardware, such as 
a convolution processor using over 2,000 correlators per 
satellite signal. Such processors have as many 

10 correlators per satellite channel as there are samples 
per code epoch, with each correlator having the same 
number of stages as the number of samples per code epoch. 
In other words, each satellite channel uses N^ correlator 
stages, where N = number of samples per code epoch. 

15 (See, e.g., the paper titled "Indoor GPS Technology" by 
Frank van Diggelen & Charles Abrams of Global Locate, 
Inc., presented at the CTIA Wireless Agenda, Dallas, May 
2001.) 

At decisional step 1434, a determination is made 
20 regarding whether a correlation peak exists in each of 
the pseudorange time sequences that is sufficiently above 
a noise -background threshold. According to one 

embodiment, a correlation peak is sufficiently above the 
noise -background threshold when the correlation peak is 
25 at least 6dB above the noise -background threshold. 
However, it will be understood that a correlation peak 
may be sufficiently above the noise-background threshold 
when the correlation peak is any suitable level above the 
noise -background threshold. If such a peak does not 
3 0 exist for any of the pseudorange time sequences, the 
method follows the No branch from decisional step 1434 to 
decisional step 1436. 
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At decisional step 1436, a determination is made 
regarding whether a maximum number of signal samples that 
have been combined for pseudorange estimation has been 
exceeded for the pseudorange time sequences without a 
5 correlation peak. According to one embodiment, the 
maximum number of signal samples that have been combined 
for pseudorange estimation comprises the number of 
samples that may be combined in approximately one second. 
Thus, the maximum number of signal samples that have been 
10 combined for pseudorange estimation may be based on the 
sampling rate. However, it will be understood that the 
maximum number of combined samples may be any suitable 
number . 

If the maximum number of combined pseudorange 

15 samples has been exceeded, this indicates that the 
corresponding pseudorange signals being received are too 
weak at the ranging receiver 14 to be used. In this 
situation, the method follows the Yes branch from 
decisional step 1436 to step 1437. At step 1437, the 

20 ranging receiver 14 sets a flag to indicate that the 
signal is too weak and returns to step 1400 where the 
ranging receiver 14 reverts to the low-power standby 
condition and waits for another trigger. However, if the 
maximum number of combined pseudorange samples has not 

25 been exceeded, the method follows the No branch from 
decisional step 1436 and returns to step 1422 to collect 
more signal samples for pseudorange estimation. 

Returning to decisional step 1434, if a correlation 
peak that is sufficiently above the noise -background 

3 0 threshold does exist for each of the pseudorange time 
sequences, the method follows the Yes branch from 
decisional step 1434 to step 1438. At step 1438, a 
multi-path correcting centroid calculation is applied to 
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determine a substantially exact location of the center of 
each of the peaks based on techniques that compensate for 
the distortion of each peak's shape by the presence of 
multi-path signals occurring near the correlation peak. 
5 At step 1442, the ranging receiver 14 sends 

pseudorange data for each of the carriers detected from 
the constellation of satellites 12 to the geolocation 
processor 24 for calculation of the geolocation of the 
ranging receiver 14. According to one embodiment, the 

10 pseudorange data comprises the estimate of the 
pseudorange, the measurement time, the satellite ID for 
each satellite 12 and/or any other suitable data. At 
step 1444, the ranging receiver 14 or other requesting 
application receives the geolocation of the ranging 

15 receiver 14 from the geolocation processor 24, after 
which the method returns to step 1400 where the ranging 
receiver 14 reverts to the low-power standby condition 
and waits for another trigger. 

FIGURE 15 is a flow diagram illustrating a method 

20 for determining geolocation in the geolocation processor 
24 by providing matching assistance in accordance with 
one embodiment of the present invention. The method 
begins at step 1500 where the geolocation processor 24 
receives geolocation initiating information. For 

25 example, the geolocation initiating information may 
comprise protocol identifying information and version 
information for the ranging receiver 14 provided by the 
ranging receiver 14 in response to a trigger. However, 
it will be understood that the geolocation initiating 

3 0 information may comprise any suitable information related 
to the initiation of the geolocation process. According 
to one embodiment, the geolocation initiating information 
is received from the ranging receiver 14 over the voice 
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channel for the wireless device 16 by way of the service 
center 18 . 

At step 1502, the geolocation processor 24 sends an 
acknowledgement of the receipt of the geolocation 
5 initiating information. For example, the geolocation 
processor 24 may transmit a handshake code to the ranging 
receiver 14 in order to establish communication and 
acknowledge compatibility. 

At step 1504, the geolocation processor 24 receives 

10 a time calibration signal from the ranging receiver 14 . 
The time calibration signal may comprise a complete 
sequence of, for example, a maximal length pseudorandom 
noise sequence or a Barker code modulated on an audio 
carrier tone. It will be understood, however, that the 

15 time calibration signal may comprise any other suitable 
timing signal. 

At step 1506, the geolocation processor 24 sends a 
time calibration response based on the time calibration 
signal. According to one embodiment, the geolocation 

20 processor 24 sends the time calibration response a pre- 
defined period of time after receiving the time 
calibration signal. For example, the pre-defined period 
of time may comprise 10 milliseconds or other suitable 
period of time. The time calibration response may 

25 comprise a different maximal length pseudorandom noise 
sequence modulated on the same carrier tone as the time 
calibration signal, followed by a message indicating a 
real time associated with the response. The real time 
information may be used by the ranging receiver 14 in 

30 order to calibrate the real-time clock 232. 

At step 1508, the geolocation processor 24 receives 
CFO estimates from the ranging receiver 14 . At step 
1510, the geolocation processor 24 matches the 
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differences in the CFO estimates from the ranging 
receiver 14 to differences in Doppler frequency shift 
measurements for signals received by the receivers 48 
and/or 6 0 or based on geometric prediction from the 
5 satellite constellation's current geometric data. 

At step 1512, any duplicate matches found in step 
1510 are eliminated in order to identify the satellites 
12. According to one embodiment, duplicate matches are 
eliminated by a least-squares regression method. At step 

10 1514, the geolocation processor 24 sends satellite 
information to the ranging receiver 14 . The satellite 
information may comprise the satellite IDs, the satellite 
message data patterns, the timing information for bit 
transitions in the satellite message data patterns and/or 

15 any other suitable information. 

At step 1516, the geolocation processor 24 receives 
pseudorange data, which is generated based on the 
satellite information, from the ranging receiver 14. The 
pseudorange data may comprise an estimate of the 

20 pseudorange, a measurement time, a satellite ID for each 
satellite 12 and/or any other suitable data. 

At step 1518, the geolocation processor 24 
determines the geolocation of the ranging receiver 14 
based on the pseudorange data, the current ephemeris data 

25 and/or any other suitable information. At step 1520, the 
geolocation processor 24 sends the geolocation of the 
ranging receiver 14 to the requesting application, at 
which point the method comes to an end. The geolocation 
may be sent in-band over the voice channel of the 

30 wireless device 16, over the data network 66, or by any 
other suitable means. 

Although the present invention has been described 
with several embodiments, various changes and 
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modifications may be suggested to one skilled in the art. 
It is intended that the present invention encompass such 
changes and modifications as fall within the scope of the 
appended claims. 
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